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EXECUTIVE  SUMMARY 


OBJECTIVES 

The  principal  objective  of  the  present  experiments  is  to  study  the  response  of  luminescent 
organisms  to  hydrodynamic  stimuli  in  a  well  characterized  flow  field.  In  particular,  using  fully 
developed  pipe  flow:  (1)  the  bioluminescent  response  of  freshly  collected  organisms  are  com¬ 
pared  in  laminar  and  turbulent  flow,  (2)  threshold  values  of  wall  shear  stress,  Twain  are  deter¬ 
mined,  and  (3)  the  flash  level  of  individual  organisms  as  a  function  of  Xwaii  are  examined. 

RESULTS 

Several  trends  describing  the  response  of  in  situ  bioluminescent  organisms  to  fully  devel¬ 
oped,  pipe  flow  stimuli  are  reported.  These  include  the  repeated  observation  that  the  biolumines¬ 
cent  organisms  studied  are  stimulated  in  laminar  flow  at  threshold  levels  of  wall  shear  stress, 
'Twain  of  the  order  of  a  few  (0.5-3)  dynes/cm^.  It  is  argued  that  this  level  of  threshold  stimulus 
should  provide  antipreditation  value.  For  the  same  value  of  Twain  it  does  not  appear  that  turbu¬ 
lent  flow  is  particularly  more  stimulatory  than  laminar  flow  on  an  individual  cell  basis  (more 
cells  may  be  excited  in  turbulent  flow  but  peak  values  do  not  appear  to  be  effected).  Another 
general  feature  of  these  measurements  is  that  with  increasing  Twaii  a  conspicuous  break  in  the 
rate  of  increase  of  average  bioluminescence  normally  occurs  between  10  to  20  dynes/cm^  (lami¬ 
nar  flow),  with  the  greater  increase  occurring  initially.  This  range  of  Twaii  is  also  where  peak 
flash  levels  of  individual  organisms  become,  after  initially  exhibiting  a  graded  response  (increas¬ 
ing  with  Twaii  for  'Twaii  <  10-20  dynes/cm^),  approximately  constant  thereafter. 

RECOMMENDATIONS 

It  is  hoped  that  this  work  will  lead  to  better  predictive  capabilities  toward  assessing  biolumi¬ 
nescent  stimulation  in  a  variety  of  oceanic  flow  fields.  Further  work  with  a  greater  range  of  Twaii 
and  Retrans  values  is  necessary  (preferably  complimented  by  laboratory  experiments  with  cul¬ 
tures  of  representative  luminous  dinoflagellates)  to  corroborate  this  initial  study.  Also  the  reac¬ 
tion  of  these  luminescent  organisms  to  varying  shear  stress  (transitional  flows)  and  acceleration 
(developing  flows)  needs  to  be  investigated. 
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1.  Photograph  of  bioluminescence.  The  marine  bioluminescent  dinoflagellate, 

Gonyaulax  polyedra  (culture  of  — 10^  cells  per  liter),  are  stimulated  as  they 

are  poured  over  a  17.5  cm  model  of  a  dolphin .  23 

2.  Computer  enhanced  photographs  of  in  situ,  flow  induced  bioluminescence 
obtained  via  a  low-light  video  camera.  Objects  can  be  identified  3-5  meters  beneath 
the  sea  surface  as  a  result  of  stimulated  bioluminescence.  At  these  depths  it  is 
likely  that  one  or  more  species  of  dinoflagellates  are  responsible  for  the 
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6b.  PMT  levels  vi'  wall  shear  stress.  San  Diego  Bay  data,  collected  67  cm  from  the 
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7a.  Representative  San  Clemente  Island  PMT  levels  vs  Reynolds  number.  Each  data 
point  represents  a  100-s  average  of  two  PMTs,  located  22  and  67  cm  from  the  inlet 
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turbulent  shear  stress,  ,  across  a  pipe  in  fully  developed  flow.  Turbulent 

shear  stress  levels  have  been  normalized  by  the  wall  shear  stress,  Twaii,  allowing 
different  Reynolds  number  experiments  to  be  compared.  The  solid  line  shows 
the  distribution  of  the  theoretical  values  of  the  total  normalized  shear  stress 

across  the  pipe.  Except  near  the  pipe  wall  (r  -+  R,  p  u  v  — »•  0;  not  shown )  the 
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Relative  PMT  level  (Imean/(Ion))  vs  Reynolds  number  as  predicted  by  equation 
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V  =  0.0104  cm^/sec  (25.6°C).  Decay  rate  of  flash  (y)  is  assumed  to  be 

0.1  sec.  Limiting  behavior  for  (Li^2/Uavg)/Y  «  1  (equation  3.11,  solid 

line)  and  (Li^2/Uavg)/Y  »  1  (equation  3.13,  dashed  line)  are  included . 


SUMMARY 


In  situ  luminous  marine  organisms  are  subjected  to  a  wide  range  of  quantifiable,  hydrody¬ 
namic  agitation  provided  by  fully-developed  laminar  and  turbulent  pipe  flow.  The  threshold  for 
stimulation  occurs  in  laminar  flow,  and  in  terms  of  the  shear  stress  at  the  pipe  wall,  Xwaih  is  only 
a  few  dynes/cm^.  With  increasing  Twaii,  a  monotonic  increase  in  the  time  averaged  biolumines¬ 
cence  is  generally  observed  throughout  the  measured  range  (Xwaii  ^  200  dynes/cm ^).  Peak  levels 
associated  with  individual  bioluminescent  flashes  behave  differently.  After  an  initial  rapid  rise 
from  threshold  to  Xwaii  values  of  10-20  dynes/cm^  (also  occurring  in  laminar  flow),  peak  levels 
remain  thereafter  nearly  independent  of  shear  stress. 

1.0  INTRODUCTION 

Awak’d  before  the  rushing  prow. 

The  mimic  fires  of  ocean  glow. 

Those  lightnings  of  the  wave; 

Wild  sparkles  crest  the  broken  tides. 

And  flashing  round,  the  vessel’s  sides 
With  elfish  lustre  lave; 

While  far  behind,  their  livid  light 
To  the  dark  billows  of  the  night 

a  blooming  splendour  gave 

from  Lord  of  the  Isles  (1815) 
by  Sir  Walter  Scott 

The  photograph^  in  figure  1  attempts  to  display  the  fairy-like  glow  of  bioluminescence.  In 
this  photograph  a  culture  of  the  bioluminescent  dinoflagellate,  Gonyaulax polyedra,  is  poured 
from  a  beaker  over  a  17.5-cm  model  of  a  dolphin.  The  concentration  of  the  organisms  was  esti¬ 
mated  to  be  of  the  order  of  10^  c/1,  representative  of  dinoflagellate  bloom  conditions  (Seliger  et 
al.,  1970).  Along  the  dolphin  model,  the  eye  and  white  underside  is  clearly  illuminated  by  the 
blue-green  light  (  —  480  nm,  Hastings  &  Sweeney  (1956);  Seliger  et  al.  (1969))  emitted  by  the 
organisms  (- 1-2  X  10^  photons/cell  (Biggley  et  al.  (1969);  Krasnow  et  al.  (1981)).  Notice  that 
some  bioluminescent  organisms  are  excited  before  striking  the  model.  It  is  not  known  whether 
the  shear  stresses  associated  with  the  pouring  motion  or  the  accompanying  acceleration  is  most 
important  towards  stimulating  bioluminescence  upstream  of  the  model.  Although  there  is  some 
evidence  (Christianson  &  Sweeney,  1972)  suggesting  the  former,  further  work  is  necessary  to 
differentiate  between  the  stimulatory  effects  of  shear  and  acceleration  on  bioluminescence. 

In  an  overview  of  coastal  bioluminescence,  Morin  (1983)  discusses  the  photic  environment 
wherein  bioluminescence  may  be  visually  observed.  At  the  sea  surface  bright  daylight  is  about 
six  orders  of  magnitude  greater  than  maximum  bioluminescent  intensities  (which  have  been 
determined  to  be  greater  than  10"^  pW  cm"^).  However,  light  levels  at  the  sea  surface  during  a 
moonless  night  may  be  two  orders  of  magnitude  less  than  maximum  bioluminescent  intensities. 

^  The  photograph  was  taken  with  Super  HG  Fujicolor  1600  film,  with  an  f  stop  of  2.8,  and  an  exposure  time  of  7 
seconds. 
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Seliger  et  al.  (1969)  have  made  the  interesting  observation  that  while  concentrations  of  one  Pyro- 
dinium  bahamense  per  milliliter  are  sufficient  to  observe  the  “white”  bioluminescent  wake  of  a 
boat,  concentrations  of  several  hundred  per  milliliter  provide  a  brilliant  glow  bright  enough  to 
read  by. 

The  computer  enhanced  photographs  of  flow-induced  bioluminescence,  shown  in  figure  2, 
were  obtained  at  sea  from  a  low-light  video  camera.  They  are  illustrative  of  the  bioluminescent 
displays  visible  to  the  naked  eye,  when  background  lighting  is  suitably  low,  even  when  the  stim¬ 
uli  is  relatively  weak  (as  compared  with  breaking  waves  or  a  ship’s  wake).  In  figures  2a  through 
c  the  configuration  of  a  rope,  the  upstream  profile  of  a  5-gallon  box,  and  the  outline  of  a  24-cm 
diameter  sphere  are  all  discernable  a  few  meters  beneath  the  sea  surface  as  a  result  of  the  biolu¬ 
minescent  light  stimulated  by  passing  swell.  The  bioluminescent  irradiance  of  the  0.25-m  diame¬ 
ter  sphere  (2c)  was  estimated  to  be  about  2.9x10^^  photons  per  second.  Figure  2d  is  representa¬ 
tive  of  the  bioluminescent  signatures  of  the  rising  bubbles  (center-left  of  photo)  and  swimming 
motion  (center-right)  of  divers  4  m  below  the  sea  surface. 

Unfortunately,  the  uncontrolled  flow  fields  associated  with  these  pictures  are  not  readily 
amenable  to  analysis  for  determination  of  the  nature  and  level  of  the  associated  hydrodynamic 
stimulus.  Moreover,  most  common  laboratory  methods  for  stimulating  bioluminescence,  e.g., 
bubbling  (Biggley  et  al.,  1969),  stirring  (Hamman  &  Seliger,  1972),  vibration  (Eckert  & 

Sibaoka,  1968),  and  chemical  (Hamman  &  Seliger,  1972)  and  electrical  (Eckert,  1966)  excita¬ 
tion,  either  do  not  provide  hydrodynamic  stimuli,  or  if  so,  can  not  provide  a  detailed  description 
thereof.  Only  through  observing  the  response  of  luminescent  organisms  in  a  well-characterized 
flow  field  (e.g.,  pipe  flow)  can  the  hydrodynamic  parameters  responsible  for  bioluminescent 
stimulation  be  best  studied.  It  is  hoped  that  this  approach  will  eventually  lead  to  reliable  esti¬ 
mates  of  the  bioluminescent  potential  of  many  of  the  flow  fields  encountered  in  the  ocean.  For 
example,  once  bioluminescence  is  determined  as  a  function  of  wall  shear  stress  from  in  situ  mea¬ 
surements,  it  should  be  possible  to  estimate  the  bioluminescence  excited  by  a  submerged  body 
moving  through  the  same  medium. 

The  objectives  of  the  present  pipe  flow  study  are  (1)  to  compare  the  response  of  biolumines¬ 
cence  in  fully  developed  laminar  and  turbulent  flow,  (2)  to  determine  threshold  values  of  wall 
shear  stress,  Twaih  necessary  to  excite  bioluminescence,  and  (3)  to  examine  the  level  of  individual 
flashes  as  a  function  of  Xwaii- 

Section  2  describes  the  experimental  apparatus  and  procedure.  Wall  shear  stress,  volume 
flow  rates,  and  photomultiplier  tube  (PMT)  measurements  are  presented  in  section  3.  These  sec¬ 
tions  focus  respectively  on:  the  response  of  luminous  organisms  before  and  after  transition  from 
laminar  to  turbulent  flow  (section  3.1a  -  San  Diego  Bay  measurements,  section  3.1b  -  San  Cle¬ 
mente  Island  measurements);  laminar  flow  stimuli  (section  3.2a  -  high  levels,  section  3.2b  - 
threshold  levels);  and  the  flash  characteristics  of  individual  organisms  (section  3.3).  A  general 
discussion  of  the  experimental  results  can  be  found  in  section  4  which  is  divided  into  subsections 
addressing:  comparison  with  previous  experiments  (section  4.1),  and  the  response  of  biolumines¬ 
cence  in  laminar  (section  4.2)  and  turbulent  (section  4.3)  flow.  Concluding  remarks  follow  in 
section  5.  Appendix  1  provides  a  general  description  of  the  methods  involved  in  estimating  the 
light  levels  associated  with  the  photographs  in  figure  2.  In  appendix  2  the  character  of  fully 
developed  laminar  and  turbulent  pipe  flow,  important  to  this  study,  is  briefly  outlined.  Finally,  in 
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appendix  3  the  mathematical  framework  is  developed  for  determining  the  bioluminescence 
dependence  on  the  convective  velocity  in  pipe  flow. 
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2.0  MATERIALS  AND  METHODS 


Fully  developed  pipe  flow  is  chosen  to  stimulate  bioluminescence  because  it  (1)  is  well  char¬ 
acterized,  (2)  offers  a  wide  range  of  laminar  and  turbulent  stimuli,  (3)  provides  continuous 
replacement  of  luminous  organisms,  and  (4)  levels  of  bioluminescence  have  not  appeared  signif¬ 
icantly  degraded  by  the  organisms  upstream  transit  through  the  pipe  (Losee  and  Lapota,  1981). 
Interest  in  pipe  flow  dates  back  to  antiquity,  with  the  earliest  pipes  being  composed  of  natural 
tubes  such  as  bamboo.  Clay  water  pipes  were  constructed  by  the  Greeks  at  least  2500  years  ago 
and  by  312  B.C.  the  water  supply  system  of  Rome  provided  as  much  water  per  capita  as  is  pres¬ 
ently  found  in  many  parts  of  the  civilized  world  (Rouse  &  Ince,  1957). 

An  exact  solution  to  the  governing  hydrodynamic  equations  for  fully  developed  laminar  pipe 
flow,  as  referenced  by  Rouse  and  Ince  (1957),  has  been  known  for  over  150  years.  Measure¬ 
ments  which  would  corroborate  these  analytical  results  had  been  already  made  (Hagen,  1839; 
Poiseuille,  1840).  In  contrast,  corresponding  knowledge  of  fully  developed  turbulent  pipe  flow  is 
predominantly  empirical  in  nature.  Nevertheless,  a  detailed  description  of  the  turbulent  flow 
field  has  emerged  (Laufer,  1954).  A  brief  development  of  the  relevant  pipe  flow  equations  used 
throughout  the  present  study  can  be  found  in  appendix  2. 

The  primary  apparatus  (figure  3)  consisted  of  a  75-liter  tank,  a  gently  contracting  inlet  sec¬ 
tion,  and  a  0.635-cm  i.d.  polycarbonate  pipe  (Rohr  et  al.,  1990).  Dye  visualization  showed  that 
the  flow  through  the  pipe  inlet  remained  laminar  throughout  the  range  of  flows  used.  Pressure 
taps  were  located  at  28.6,  60.4,  and  92.2  cm  (or  about  45,  95,  and  145  pipe  diameters)  from  the 
pipe  inlet.  Pressure  drops  were  measured  by  variable  reluctance  differential  transducers  (Valy- 
dyne  Corporation).  No  significant  difference  was  found  between  pressure-drop  measurements 
across  the  first  and  second  taps,  and  the  second  and  third  taps  (supporting  the  premise  that  the 
flow  is  fully  developed).  Wall  shear  stress  values  were  calculated  from  pressure  drop  measure¬ 
ments  employing  the  second  and  third  pressure  taps  (see  equation  2.1b). 

On  one  occasion,  for  the  purpose  of  achieving  higher  levels  of  laminar  shear  stress,  a  similar 
pipe-flow  apparatus,  but  with  a  0.305-cm  i.d.  glass  pipe,  was  employed.  The  entrance  section 
consisted  of  a  conical  section  that  tapered  from  a  5.08-cm  i.d.  to  a  0.305-cm  i.d.  over  a  length  of 
7.62  cm.  Pressure  measurement  locations  were  2.41  m  and  3.02  m  from  the  pipe  entrance.  Here 
glass  was  extruded  from  the  pipe,  cut,  and  fused  to  a  graded  glass  seal  containing  a  metal  pres¬ 
sure  tap. 

Volume  flow  rates  were  measured  at  the  discharge  end  of  each  pipe  using  either  a  rising-ball 
flow  meter  or  a  mass  flow  meter  (Micro  Motion  Corporation  Model  D40).  To  increase  accuracy 
at  low  flow  rates,  the  volume  flux  was  also  measured  by  weighing  the  volume  of  water  collected 
at  the  pipe  outlet  over  a  measured  time.  The  average  flow  velocity,  Uavg,  was  obtained  by  divid¬ 
ing  the  volume  flow  rate  by  the  pipe  area.  Measurements  of  volume  flow  rate,  pipe  diameter,  and 
temperature  (from  which  the  kinematic  viscosity  can  be  determined)  were  used  to  determine  the 
corresponding  Reynolds  number  (see  equation  2.3). 

Figure  4  shows  representative  Darcy  friction  factor  (see  equation  2.5)  vs  Reynolds  number 
data.  For  the  0.635-cm  i.d.  pipe  and  Re  <  2000,  there  is  excellent  agreement  between  the  theoret¬ 
ical  (solid  line)  and  measured  values.  Between  Re  >  2000  and  transition  from  laminar  to 
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turbulent  flow,  Retrans^  the  friction  factor  data  for  the  0.635-cm  i.d.  pipe  is  somewhat  higher  than 
the  theoretical  predictions.  This  discrepancy  results  because  the  pressure  taps  for  the  0.635-cm 
pipe  were,  in  terms  of  pipe  diameter,  not  far  enough  downstream  (Lg/D  =  95  and  145,  see  equa¬ 
tion  2.2a)  to  ensure  fully  developed  laminar  flow  for  this  range  of  Reynolds  numbers.  Corre¬ 
sponding  Twaii  calculations  will  be  higher,  about  40%  at  Re  =  6000,  than  what  would  be  pre¬ 
dicted  for  fully  developed  laminar  flow.  The  pressure  taps  for  the  0.305-cm  i.d.  pipe,  however, 
were  located  much  further  downstream  (LJD  =  960).  Consequently,  better  agreement  between 
measurements  and  predicted  values  are  found  for  the  0.305-cm  pipe  at  the  higher  Reynolds  num¬ 
ber  laminar  flows.  Figure  4  also  shows,  for  Re  >  Retrans^  excellent  agreement  between  the 
accepted  turbulent  values  (Blasius,  1913;  dashed  line)  and  the  present  measurements.  Simulta¬ 
neous  friction  factor  and  Reynolds  number  values  were  calculated  for  each  bioluminescent  mea¬ 
surement  throughout  this  investigation.  This  procedure  allowed  both  the  nature  of  the  flow  field 
(i.e.,  whether  it  was  laminar  or  turbulent)  to  be  determined,  and  shear  stress  levels  throughout  the 
pipe  to  be  quantified. 

Bioluminescence  was  measured  by  RCA  8575  photomultiplier  tubes  (PMTs)  used  in  the  pho¬ 
ton  count  mode.  The  face  of  each  PMT  fits  into  an  opaque  housing  which  clamped  directly  onto 
the  pipe.  With  this  arrangement,  the  PMT  field  of  view  encompassed  5.0  cm  of  pipe  length.  For 
the  0.635-cm  i.d.  pipe,  the  center  of  one  PMT  was  always  positioned  near  the  bottom  of  the  pipe, 
67  cm  («*  106  pipe  diameters)  from  the  inlet  (figure  2).  A  second  PMT  was  sometimes 
employed,  positioned  either  near  the  top  of  the  pipe  with  its  center  6  cm  ( «  9  pipe  diameters) 
from  the  inlet,  or  near  the  middle  of  the  pipe  with  its  center  22  cm  ( *=  35  pipe  diameters)  from 
the  inlet.  For  the  0.305-cm  i.d.  glass  pipe  the  the  PMT  center  was  located  1.06  m  (~  350  pipe 
diameters)  from  the  inlet. 

Throughout  these  experiments,  surface  water  from  San  Diego  Bay  and  off  San  Clemente 
Island  was  obtained  by  bucket  after  dusk.  On  two  occasions,  luminescent  organisms  that  passed 
through  the  pipe  were  identified.  The  majority  of  organisms  found  were  Gonyaulax polyedra  and 
Noctiluca  millepias.  These  dinoflagellates  are  common  to  the  coast  of  Southern  California 
(Sweeney,  1963;  Holmes  et  al.,  1967;  Kimor,  1983).  Together  with  the  genera  Ceratium  and  Pro- 
toperidinium  they  are  expected,  to  various  degrees  throughout  the  year,  to  be  the  principal 
sources  of  shallow-water  bioluminescent  light  occurring  in  San  Diego  Bay  (Sweeney,  1963; 
Kimor,  1983). 

Most  bioluminescent  dinoflagellate  flashes  last  about  100  ms  (Anderson  et  al.,  1988).  For 
this  flash  duration,  streaks  of  bioluminescent  light  will  be  less  than  25  cm  long  at  the  highest 
flow  rates  achieved  in  the  present  experiments.  Therefore,  the  light  contribution  to  the  bottom 
PMT  from  flashes  initiated  at  the  pipe  inlet,  67  cm  upstream,  are  presumed  inconsequential. 

Even  allowing  for  a  100  ms  response  time  and  a  200  ms  flash  duration,  as  shown  in  figure  5, 
light  streaks  initiated  at  the  inlet  are  not  measured  by  the  bottom  PMT  throughout  the  flow  range 
of  interest.  The  lack  of  significant  influence  of  inlet  stimulus  on  the  bottom  (67  cm  from  inlet) 
PMT  measurements  is  corroborated  by  recordings  from  a  low-light  video  camera  (Cohu  Inc. 
ISIT)  with  a  video  analyzer  (Hamamatsu  model  DVS-3000).  Finally,  lack  of  time  correlation 
between  the  top  and  bottom  PMTs  also  indicates  that  light  is  not  being  internally  reflected  down 
the  pipe. 
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The  average  PMT  measurements  presented  here  are  in  the  form  of  counts  per  second  (Morita 
&  Johnson,  1974;  Rohr  et  al.,  1990),  as  opposed  to  counts  per  volume  (Widder  et  al.,  1993;  Latz 
et  al.,  1994).  Threshold  levels  of  Xwaii  are  better  displayed  when  PMT  levels  are  expressed  as 
counts  per  second.  It  is  also  more  appropriate  to  express  individual  flash  PMT  levels  as  counts 
per  second  as  opposed  to  counts  per  volume.  However,  expressing  PMT  levels  as  counts  per 
second  may  obscure  the  interpretation  of  why  average  PMT  levels  change  with  increasing  flow 
rate.  An  inherent  ambiguity  results  because  both  convective  and  stimulative  effects  determine 
PMT  levels.  In  appendix  3  the  dependence  of  the  average  PMT  count  rate  on  the  convective 
velocity,  for  the  present  pipe  flow  configurations,  is  shown  to  be  approximately  linear  through¬ 
out  the  flow  regime  examined.  Therefore,  accompanying  each  PMT  vx  Xwaii  curve  presented, 
where  the  PMT  levels  are  expressed  as  counts  per  second,  the  corresponding  relation  between 
PMT  levels  expressed  as  counts  per  volume  (calculated  via  equations  2.6b  and  2.11)  and  Xwaii  is 
also  included. 

The  entire  pipe-flow  apparatus  resided  within  a  dark  room  and  background  PMT  levels  were 
typically  less  than  200  counts  per  second.  Average  pressure  drop,  PMT,  and  flow  rate  measure¬ 
ments  were  obtained  for  either  equal  flow  volumes  (2.2  liters  for  10/1, 10/91  experiments)  or 
equal  times  (100  seconds,  remaining  experiments),  the  results  being  essentially  the  same.  The 
time  series  for  the  PMT  data,  containing  1000  consecutive  0.1 -s  increments,  were  also  recorded. 
The  time  series  were  used  as  a  diagnostic  to  discriminate  between  background  noise  and  actual 
bioluminescence,  and  to  provide  information  about  the  level  of  individual  bioluminescent 
flashes. 
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3.0  RESULTS 


3.1  BIOLUMINESCENCE  BEFORE  (LAMINAR)  AND  AFTER  (TURBULENT) 
TRANSITION 

The  following  experiments,  which  employed  the  0.635-cm  i.d.  pipe  (see  figure  3),  focus  on 
the  response  of  bioluminescence  on  either  side  of  transition  from  fully  developed  laminar  to 
fully  developed  turbulent  pipe  flow.  During  transition,  the  flow  was  characterized  by  the  inter¬ 
mittent  passage  of  laminar  and  turbulent  slugs  of  water,  which  resulted  in  large  fluctuations  in 
the  velocity  and  pressure  measurements.  Transitional  flow  measurements  are  for  the  most  part 
(except  some  San  Clemente  Island  measurements  -  see  figure  7c)  not  included  in  the  present 
study. 

3.1a  San  Diego  Bay  Measurements 

Bioluminescent  organisms  collected  from  San  Diego  Bay  surface  waters  were  allowed  to 
dark  adapt  in  the  tank  for  at  least  45  minutes  before  beginning  the  pipe  flow  experiments.  Spe¬ 
cific  times  marking  the  commencement  of  data  collection  are  included  throughout  the  text  and 
figure  legends.  Laminar  to  turbulent  transition  occurred  between  Reynolds  numbers  of  5000  to 
6000.  The  position  of  the  single  PMT  for  these  measurements  was  67  cm  (106  pipe  diameters) 
from  the  pipe  entrance. 

Representative  San  Diego  Bay  PMT  measurements  are  presented  in  figure  6a  as  a  function  of 
Reynolds  number.  Each  experiment  involved  a  single  tank  ( —  65  liters  full)  of  bay  water  and  was 
divided  into  two  runs.  The  open  symbols  in  figure  6a  represent  measurements  from  the  initial 
run,  designated  run  1,  which  employed  about  half  of  the  reservoir.  The  sequence  of  flow  rates 
was  monotonically  decreasing  during  run  1.  The  solid  symbols  in  figure  6a  are  measurements 
collected  from  the  final  run,  designated  run  2,  which  comprised  the  remaining  half  of  the  reser¬ 
voir.  Here  the  sequence  of  flow  rates  was  continually  increasing.  Each  data  point  in  figure  4a 
represents  an  average  over  2.2  liters  of  seawater. 

An  increasing  difference  between  PMT  values  associated  with  the  first  and  second  runs,  for 
nearly  the  same  Reynolds  number,  is  observed  in  figure  6a  for  both  experiments.  This  divergence 
between  runs  1  and  2  may  reflect,  due  to  the  different  tank  residency  times  between  runs,  the  dif¬ 
ferent  prestimulus  recovery  times  of  the  luminous  organisms.  Biggley  et  al.  (1969)  have  found 
Gonyaulax  polyedra,  during  their  dark  cycle,  to  recover  exponentially  after  being  stimulated  to 
exhaustion  with  half  times  of  about  20  minutes.  While  the  stimulus  accompanying  the  collection, 
transportation,  and  tank  filling  is  presumably  much  less  than  that  associated  with  total  depletion, 
a  similar  “hysteresis”  effect,  as  shown  in  figure  6a,  is  normally  evident  throughout  the  data 
where  this  flow-rate  sequence  is  followed.  Another  possible  explanation  for  this  apparent  tank 
residency  time  dependence  is  that  stratification  of  the  luminous  organisms  had  either  occurred 
prior  or  during  the  experiment.  Also  apparent  in  figure  6a  is  an  increase  in  the  slope  of  the  PMT 
vs  Re  data  around  transition  (5000  <  Retrans  —  6000).  This  behavior  was  normally  observed 
and  often  more  dramatic  for  larger  values  of  Rctrans- 

Figure  6b  plots  the  same  PMT  measurements  of  figure  6a,  but  as  a  function  of  the  corre¬ 
sponding  wall  shear  stress,  Xwaii-  Included  in  figure  6b  are  additional  measurements  taken  from 
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San  Diego  Bay  throughout  the  same  flow  regime  and  with  similar  Retrans  values.  (The  cross 
symbols  are  comparison  data  from  Rohr  et  al.  (1990)  which  are  addressed  in  section  4.)  The 
abrupt  increase  in  bioluminescence  appearing  at  transition  in  figure  6a  is  no  longer  evident  in 
figure  6b. 

The  solid  line  in  figure  6b  represents  the  slope  (1.03,  R^  =  0.49)  of  the  least-squares  power 
law  regression  for  all  the  San  Diego  Bay  data  for  Xwaii  >  25  dynes/cm^  (turbulent  flow).  Dividing 
the  PMT  measurements  in  figure  6a  by  the  corresponding  volume  flux  (via  eq.  2.11)  results  in  a 
slope  of  0.46,  i.e.,  PMT  levels  expressed  as  counts  per  volume  oc  Vaii®'^^.  As  shown  in  appen¬ 
dix  3,  this  latter  representation  eliminates  the  flow’s  convective  contribution  towards  increasing 
PMT  levels. 

3.1b  At-Sea  Measurements  Off  San  Clemente  Island 

To  complement  the  San  Diego  Bay  studies,  the  same  pipe-flow  apparatus  was  used  in  a  series 
of  at-sea  experiments  conducted  from  the  research  vessel  USNS  De  Steiguer.  The  experiments 
were  performed  a  few  miles  off  the  southeast  coast  of  San  Clemente  Island  (32°  47.6-56.1'  N, 
118°  25.1 '-25.5'  W)  from  3/7/92  to  3/11/92.  In  an  attempt  to  mitigate  the  bias  introduced  by  dif¬ 
ferent  reservoir  residency  times  of  the  organisms,  dark  adaption  times  were  reduced  to  20  min¬ 
utes  and  the  at-sea  flow-rate  sequence  alternated  between  high  and  low  flow  rates. 

Two  PMTs  were  employed  throughout  these  at-sea  experiments.  The  PMTs  were  located  22 
cm  (middle  position)  and  67  cm  (bottom  position)  from  the  pipe  inlet  (see  figure  3).  Measure¬ 
ments  from  the  two  PMTs  usually  exhibited  good  agreement  and  consequently  were  averaged. 
Transitional  Reynolds  numbers  were  found  to  vary  significantly  throughout  this  series  of  experi¬ 
ments  (Retrans  ~  3500,  5400,  and  7600  in  figure  7a).  Presumably,  this  was  due  to  the  changing 
background  motion  accompanying  the  evolving  sea  state.  The  most  conspicuous  differences 
between  the  PMT  vs  Re  data  sets  shown  in  figure  7a  occurred  between  Reynolds  numbers  of 
3000  and  9000.  For  example,  throughout  this  range  the  PMT  levels  associated  with  low  Retrans 
were  significantly  higher  than  those  corresponding  to  high  Retrans-  This  disparity  reflects  the  two 
different  states  of  the  flow,  i.e.,  turbulent  (higher  PMT  values)  and  laminar,  existing  for  the  same 
Re  number  because  of  the  different  Retrans  values.  Similar  arguments  can  be  made  for  the  behav¬ 
ior  of  the  middle-of-the-range  Retrans  data  (i.e.  Retrans  ~  5400).  As  before,  when  comparing 
laminar  and  turbulent  flows  at  similar  Reynolds  numbers  in  figure  7a,  the  PMT  values  associated 
with  turbulent  flow  are  usually  greater. 

Plotting  the  same  average  PMT  values  in  figure  7a  as  a  function  of  wall  shear  stress  in  figure 
7b,  both  laminar  and  turbulent  flow  data  appear  to  lie  along  a  single  curve.  Note  that  the  highest 
laminar  Xwaii  value  (3/9/92,  A;  Re  =  7600,  Xwaii  ~  35  dynes/cm^)  overlaps  the  lowest  turbulent 
Xwaii  value  (3/9/91-3/10/92,  □;  Re  »  4500,  Xwaii  ~  30  dynes/cm^).  Also  noteworthy  is  that 
within  the  laminar  flow  region  of  each  data  set  in  figure  7b,  a  change  of  slope  around  Xwaii  ~ 
15-20  dynes/cm^  is  apparent.  This  feature,  also  evident  in  figure  6b,  will  be  addressed  in  the  dis¬ 
cussion  section.  For  x^aii  >  25  dynes/cm^,  the  San  Clemente  Island  PMT  laminar  and  turbulent 
data  in  figure  7b  increase  approximately  as  the  first  power  (slope  =  0.82,  R^  =  0.72)  of  wall 
shear  stress.  (Expressing  these  PMT  levels  as  counts  per  volume  results  in  a  slope  of  0.25.)  The 
solid  line  appearing  in  figure  7b,  included  for  comparison  purposes,  has  a  slope  of  1.03  as  deter¬ 
mined  from  the  San  Diego  Bay  data  (figure  6b). 
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In  figure  7c  the  previous  data  from  figure  7b  are  replotted  to  show  the  overlap  between 
experiments  of  the  different  flow  regimes.  For  all  data  sets  laminar  flow  is  indicated  by  the  open 
square  symbol,  turbulent  flow  data  by  the  solid  square  symbol,  and  transitional  flow  by  an  open 
box  with  X  inside.  In  general  the  response  of  bioluminescence  to  wall  shear  stress  in  pipe  flow,  as 
shown  in  figure  7c,  does  not  appear  to  be  significantly  dependent  on  the  nature  of  the  flow.  How¬ 
ever,  two  caveats  should  be  mentioned.  First,  unlike  the  transitional  data  appearing  in  figure  7c, 
there  have  been  similarly  acquired  measurements  collected  in  San  Diego  Bay,  which  show  a  sig¬ 
nificant  increase  in  bioluminescence  during  transition.  Secondly,  although  equation  2.1b  was 
invoked,  neither  the  transitional  nor  the  higher  laminar  flow  (Re  >  2000)  measurements  were 
obtained  in  fully  developed  flow. 

3.2  BIOLUMINESCENCE  IN  LAMINAR  FLOW 

A  rise  in  PMT  levels  with  increasing  laminar  flow  is  evident  to  various  degrees  throughout 
figures  6  and  7.  The  present  section  focuses  on  the  effect  on  bioluminescence  of  both  high  and 
low  Twaii  levels  of  laminar  flow  stimuli.  High  Xwaii  levels  obtained  in  laminar  flows  may  overlap 
low  Twaii  values  obtained  in  turbulent  flows,  thereby  providing  an  opportunity  to  compare  the 
stimulative  effect  of  laminar  and  turbulent  flows  for  the  same  average  shear  field.  Low  levels  of 
Twali  in  laminar  flow  are  important  toward  determining  threshold  values  of  the  hydrodynamic 
stimulation  of  bioluminescence. 

3.2a  High  Levels  of  Laminar  Flow  Stimulus,  Xwaii  >  15  dynes/cm^ 

The  highest  laminar  flow  Reynolds  number  (Rctrans  ~  10,000)  measurements  were  obtained 
during  the  USNS  De  Steiguer  experiments  conducted  off  San  Clemente  Island  (32°  55.1'  N, 

118°  27.7'  W)  on  12/18/91  and  12/19/91.  The  0.635-cm  i.d.  pipe-flow  apparatus  was  again  uti¬ 
lized  with  the  PMT  positioned  67  cm  from  the  pipe  inlet.  No  significant  differences  were 
observed  between  the  six  consecutive  runs  executed.  Consequently,  all  PMT  vs  Xwaii  data  were 
plotted  collectively  in  figure  8. 

The  highest  laminar  flow  shear  stress  data  were  obtained  in  the  0.305-cm  i.d.  pipe-flow  appa¬ 
ratus  with  seawater  collected  from  San  Diego  Bay  on  12/12/91.  The  PMT  was  106  cm  from  the 
pipe  inlet.  The  smaller  diameter  pipe,  while  providing  a  greater  range  of  fully  developed  laminar 
flow  than  the  0.635-cm  pipe,  significantly  reduced  the  volume  flux  for  equal  x^aii  values.  This 
resulted  in  lower  signal  levels  and  greater  variability.  To  reduce  variability,  each  data  point  from 
the  0.305-cm  i.d.  pipe  in  figure  8  is  composed  of  an  average  of  five  consecutive  100-s  records. 
For  comparison  with  the  0.635-cm  i.d.  data  in  figure  8,  the  PMT  levels  obtained  from  the 
0.305-cm  i.d.  glass  pipe  are  adjusted  by  the  ratio  of  their  areas.  No  discernable  time  dependence 
was  found  on  the  sequence  in  which  the  five  consecutive,  constant  flow-rate  measurements  were 
taken.  However,  this  was  not  the  case  for  data  sets  obtained  when  tank  residency  times  were 
much  different.  For  example,  the  last  data  set  taken,  which  had  the  highest  tank  residency  time, 
was  conspicuously  high  (Xwaii  =  45  dynes/cm^,  PMT  =  44,000  counts  per  second)  and  conse¬ 
quently  is  not  included  in  figure  8. 

Representative  laminar  and  turbulent  flow  data  from  figure  6b  (10/1/91)  are  also  shown  in 
figure  8.  There  is  a  conspicuous  gap  between  Xwaii  ~  12  and  40  dynes/cm^  for  the  10/1/91  data, 
where  transition  occurs,  which  is  spanned  by  the  high  shear  laminar  flow  data.  The  lines  in 
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figure  8  represent  the  slopes  of  a  least-squares  power  law  regression  to  each  data  set  where  Twaii 
>  15  dynes/cm^.  These  slopes  are  1.00  (laminar  flow;  =  0.71;  12/18/91-12/19/91)  for  the 
large  dashed  line,  1.26  (laminar  flow;  =  0.89;  12/12/91)  for  the  solid  line,  and  1.10  (turbulent 
flow;  R2  =  0.96;  10/01/91)  for  the  small  dashed  line.  Expressing  the  PMT  levels  as  counts  per 
volume  results  in  respective  slopes  of  0.00,  0.26,  and  0.53. 

3.2b  Low  Levels  of  Laminar  Flow  Stimulus,  x  <  15  dynes/cm^ 

Figure  9  displays  PMT  levels  as  a  function  of  wall  shear  stress  from  the  most  extensive  low 
Reynolds  number  experiments.  These  experiments  used  San  Diego  Bay  water  in  the  0.635-cm 
pipe  with  the  PMT  67  cm  downstream  from  the  pipe  inlet.  Threshold  values  of  wall  shear  stress 
in  figure  7  occur  around  1  dyne/cm^.  Corresponding  Reynolds  numbers  are  less  than  1000.  At 
these  low  Reynolds  numbers  the  flow  is  unquestionably  laminar.  Similar  experiments  (not 
shown)  throughout  the  year  corroborate  that  threshold  values  always  occur  in  laminar  flow  and 
at  levels  usually  between  0.5  and  3  dynes/cm^. 

The  increase  of  bioluminescence  with  Xwaii  shown  in  figure  9  is  much  faster  than  what  is 
found  (figures  6b,  7b,  and  8)  for  Xwaii  >  15  dynes/cm^,  where  the  relationship  is  approximately 
linear.  The  solid  line  in  figure  9  represents  a  least  squares  power  law  regression  to  the  PMT  data 
between  0.7  and  14.0  dynes/cm^  and  has  a  slope  of  2.57,  R^  =  0.59.  (Accounting  for  the  convec¬ 
tive  velocity  component,  see  appendix  2,  the  bioluminescent  levels  expressed  per  unit  volume 
are  proportional  to  Xwaii^'^^-)  The  corresponding  Reynolds  numbers  for  these  data  are  always  less 
than  5000,  and  since  Retran.s  for  these  experiments  occurs  around  6000,  all  the  data  in  figure  9 
were  obtained  in  laminar  flow. 

The  associated  f  vs  Re  measurements  for  much  of  the  data  found  in  figures  8  and  9  can  be 
found  in  figure  4.  The  corresponding  PMT  vs  Re  data  are  not  presented  because,  in  fully  devel¬ 
oped  laminar  flow,  Xwaii  is  proportional  to  Re. 

3.3  INDIVIDUAL  FLASH  LEVELS 

Representative  PMT  time  series  data  are  shown  in  figure  10a  for  Xwaii  =  15  dynes/cm^  (Re  = 
5000,  laminar  flow)  and  figure  10b  for  Xwaii  =  112  dynes/cm^  (Re  =  9950,  turbulent  flow).  The 
PMT  is  located  67  cm  from  the  inlet  of  the  0.635-cm  i.d.  pipe.  Measurements  are  composed  of 
PMT  counts  in  consecutive  0.1-s  bins.  Each  vertical  spike  in  figures  10a  and  b  is  believed  to  be 
indicative  of  a  single  flashing  organism.  Considering  that  the  corresponding  100-s  average  PMT 
level  for  figure  10a  is  about  five -fold  larger  than  for  figure  10b,  it  is  remarkable  that  the  peak 
0.1-s  PMT  level  for  either  figure  remains  nearly  the  same  (within  10%). 

Representative  peak  levels  from  the  bottom  PMT  (located  67  cm  from  the  inlet  of  the 
0.635-cm  i.d.  pipe)  are  plotted  in  figure  11  as  a  function  of  Xwaii  for  both  San  Clemente  Island 
and  San  Diego  Bay  experiments.  Peak  PMT  levels  are  obtained  by  taking  the  highest  count  in 
any  0.1-s  bin  within  each  100-s  record  (constant  Xwaii)-  To  compare  with  average  PMT  values 
which  are  in  counts  per  second,  peak  values  in  figures  11  and  12  are  divided  by  0.10. 

Transition  from  laminar  to  turbulent  flow  for  all  the  experiments  in  figure  11  was  between 
Reynolds  numbers  of  5000  and  6000.  This  range  of  Rctrans  resulted  in  the  gap  in  the  shear  stress 
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data  between  20  to  35  dynes/cm^.  Between  threshold  and  Xwaii  values  of  about  10  dynes/cm^, 
where  the  flow  is  always  laminar,  peak  PMT  levels  grow  rapidly  with  increasing  shear  suggest¬ 
ing  a  graded  bioluminescent  response.  For  Xwaii  >  10  dynes/cm^,  which  includes  both  the  highest 
levels  of  laminar  shear  stress  and  the  entire  range  of  turbulent  shear  stress,  peak  PMT  levels 
remain  nearly  constant^.  Included  in  figure  11  are  representative  100-s  average  bottom  PMT  lev¬ 
els.  An  abrupt  change  in  functional  dependence  of  both  peak  (0.1  s)  and  average  (100  s)  PMT 
levels  around  Xwaii  “  10-20  dynes/cm^  is  evident. 

Peak  top  (6  cm  from  pipe  inlet)  and  bottom  (67  cm  from  pipe  inlet)  PMT  data  sets  taken 
from  the  0.635-cm  i.d.  pipe  are  compared  in  figure  12  as  a  function  of  Reynolds  number.  With 
increasing  Reynolds  number  both  top  and  bottom  peak  levels  initially  grow,  followed  by  a  flow 
regime  where  they  assume  a  nearly  constant  value.  The  onset  of  bioluminescence  for  the  top 
PMT  is  seen  to  occur  at  a  lower  Reynolds  number,  which  is  in  agreement  with  the  100-s  average 
PMT  measurements  and  visual  observations. 

The  0.1-s  peak  PMT  data  in  figure  12  are  not  plotted  as  a  function  of  Xwaii,  because  Xwaii  is 
not  known  at  the  location  of  the  top  PMT  where  the  flow  is  not  fully  developed.  Nevertheless, 
for  the  same  Reynolds  number  gross  comparisons  between  the  flow  fields  opposite  the  top  and 
bottom  PMTs  can  be  made.  For  example,  via  dye  visualizations  the  flow  past  the  top  PMT  was 
found  to  remain  laminar  throughout  the  entire  range  of  Reynolds  numbers  investigated,  whereas 
R^trans  *=  5000-6000  for  the  bottom  PMT.  Comparing  top  and  bottom  peak  PMT  data  for  Re  > 
6000,  it  appears  that  the  nature  of  the  flow  field  (laminar  vs  turbulent)  is  not  particularly  impor¬ 
tant  towards  determining  flash  levels  of  individual  bioluminescent  organisms. 


2  This  constant  value  is  about  half  the  level  where  the  PMT  signal  becomes  clipped.  Similar  behavior,  i.e.,  peak 
PMT  level  independence  at  twaii  values  from  about  10  -  200  dynes/cm^,  is  found  when  the  experiment  is  repeated 
with  a  10“^  filter. 
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4.0  DISCUSSION 


4.1  COMPARISON  WITH  PREVIOUS  EXPERIMENTS 

There  have  been  several  previous  bioluminescent  experiments  that  are  relevant  to  the  present 
study.  While  these  experiments  did  not  usually  include  direct  measurements  of  shear  stress,  for 
comparison  sake,  estimates  of  twaii  are  made  when  possible. 

4.1a  Christianson  and  Sweeney  (1972) 

Christianson  and  Sweeney  (1972)  have  studied  bioluminescent  stimulation  of  Gonyaulax 
polyedra  cultures  in  a  0.86-mm  i.d.  coiled  capillary  tube  through  a  range  of  laminar  flows.  A 
PMT  monitored  the  emitted  light  from  the  entire  coil,  including  its  centrally  located  inlet.  They 
noted  that  luminescence  was  observed  to  take  place  principally  within  the  coil  as  the  cells  passed 
through  the  capillary.  Assuming  that  the  shear  stress  in  the  coiled  tube  is  approximately  equal  to 
that  in  a  straight  tube  with  an  equivalent  diameter^,  Twaii  values  can  be  estimated  (equation  2.6b) 
and  are  included  parenthetically.  Christianson  and  Sweeney  (1972)  observed  night  phase 
luminescence  (  a  counts  per  second)  to  initially  increase  faster  than  the  flow  rate,  but  then  from 
about  7  ml/min  (20  dynes/cm^)  to  12  ml/min  (35  dynes/cm^),  the  increase  was  nearly  propor¬ 
tional  to  the  flow  rate.  A  triggering  flow  rate  of  0.8  ml/min  (2.3  dynes/cm^)  was  inferred  by 
extrapolating  to  the  lower  flow  rate  data. 

The  above  results  of  bioluminescent  excitation  are  in  good  agreement  with  those  reported 
here.  Moreover,  Christiansen  and  Sweeney  (1972)  hypothesized  that  the  change  in  the  rate  at 
which  the  average  bioluminescence  grows  with  shear  stress  may  be  related  to  an  initial  gradation 
of  peak  flash  levels.  This  hypothesis  is  supported  by  the  present  data. 

4.1b  Anderson,  Nosenchuck,  Reynolds  and  Walton  (1988) 

Anderson  et  al.  (1988)  repeated  the  experiments  of  Christianson  and  Sweeney  (1972),  using 
the  same  coiled  capillary  tube  and  cultures  of  Gonyaulax  polyedra.  An  image  intensifier  was 
employed  to  image  the  entire  capillary  coil.  Anderson  et  al.  (1988)  found,  contrary  to  Christian¬ 
son  and  Sweeney’s  (1972)  visual  observations,  that  bioluminescence  was  not  stimulated  in  the 
uniform  section  of  the  capillary  tube,  but  only  at  the  abrupt  transition  at  its  entrance.  This  result, 
together  with  additional  measurements  obtained  in  a  converging-diverging  capillary  tube,  led 
Anderson  et  al.  (1988)  to  conclude  that  bioluminescence  is  stimulated  principally  by  changes  in 
shear,  acceleration,  and  pressure. 

The  present  experiments  also  indicate  greater  (2-3  times  downstream  values)  bioluminescent 
stimulus  occurring  near  the  pipe  inlet.  However,  stimulation  is  now  found  to  continue  throughout 

^  Strictly  speaking  the  friction  factor  in  a  curved  pipe  is  often  greater  than  that  of  an  equivalent  straight  pipe  at 
identical  Reynolds  number.  In  a  coiled  tube  the  presence  of  centrifugal  forces  result  in  a  secondary  (i.e.,  radial) 
flow  pattern  which  is  a  function  of  the  Dean  number,  De,  defined  as  follows: 

De  =  Re*(D/d)l/2 

Here  d  is  the  diameter  of  the  pipe  coil,  and  as  before,  Re  is  the  pipe  Reynolds  number,  and  D  is  the  pipe  diameter. 
For  Dean  numbers  less  than  10,  the  secondary  flow  is  small  and  the  friction  factor  is  nearly  equivalent  to  the 
corresponding  straight  pipe  values.  For  larger  Dean  numbers  the  friction  factor  of  the  coil  exceeds  that  of  an 
equivalent  straight  pipe.  The  Dean  numbers  in  the  Christian  and  Sweeney  (1972)  experiments  ranged  from  0  to 
about  35  (Uavg  =  0-12  ml/sec). 
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the  pipe  when  threshold  shear  stress  levels  are  exceeded.  Note  in  particular,  contrary  to  Ander¬ 
son  et  al.’s  (1988)  findings,  bioluminescent  stimulation  is  also  observed  in  laminar  flow  even 
though  the  shear  stress  is  constant  along  the  organism’s  trajectory.  The  reason  for  this  discrep¬ 
ancy  is  not  known. 

4.1c  Morita  and  Johnson  (1974) 

Morita  and  Johnson  (1974)  have  performed  pipe  flow  studies  in  a  0.71-cm  i.d.  glass  pipe, 
also  with  cultures  of  Gonyaulax polyedra.  Transition  was  assumed  to  occur  between  a  Reynolds 
number  of  2000  and  4000.  Estimates  of  twaii  were  made  assuming  laminar  flow  for  Re  <  2000 
(equation  2.6b)  and  turbulent  flow  for  Re  >  4000  (equation  2.9b).  Morita  and  Johnson  (1974) 
found  that  for  T^aii  >  25  dynes/cm^,  which  was  assumed  to  be  fully  developed  turbulent  flow,  the 
stimulated  emission  ( a  counts  per  second)  appeared  to  be  linearly  related  to  shear.  This  obser¬ 
vation  is  in  good  agreement  with  the  present  measurements  (figures  4b  and  5b).  However,  at 
lower  shear  stress  levels,  their  results  are  contradictory  to  those  reported  here.  Morita  and  John¬ 
son  (1974)  reported  more  light  at  twaii  =  2-3  dynes/cm^  than  at  Xwaii  =  10-15  dynes/cm^.  It  was 
hypothesized  that  this  anomalous  behavior  may  have  resulted  because  the  flash  signature  (i.e., 
the  rise  time  and/or  decay  rate  of  the  stimulable  light  emission)  is  dependent  on  stimulation 
intensity. 

4.1d  Latz,  Case  and  Gran  (1994) 

Most  recently  Latz  et  al.  (1994)  reported  a  study  of  the  excitation  of  cultures  of  biolumines¬ 
cent  dinoflagellates  (including  Gonyaulax  polyedra)  by  laminar  and  turbulent  Couette  flow.  The 
general  PMT  vi'  x  behavior  they  observed,  for  shear  stresses  between  0.3  and  16  dynes/cm^,  is  in 
excellent  agreement  with  the  present  pipe  flow  measurements.  In  particular,  they  have  found 
laminar  flow  excitation  of  bioluminescence,  threshold  shear  values  of  about  1  dyne/cm^,  a 
graded  response  of  individual  flashes  at  low  shear  stress  levels,  and  a  lack  of  sensitivity  to  the 
nature  (laminar  or  turbulent)  of  the  hydrodynamic  stimuli. 

4.1e  Rohr,  Losee  and  Hoyt  (1990) 

The  at-sea  (Sea  of  Cortez)  bioluminescent  experiments  of  Rohr  et  al.  (1990)  employed  essen¬ 
tially  the  same  0.635-cm  i.d.  pipe  apparatus  (with  the  exception  of  pressure  taps)  and  procedure 
as  used  here.  Rohr  et  al.  (1990)  reported: 

1.  an  often  abrupt  increase  in  bioluminescence  as  a  function  of  Re  number  around 
transition, 

2.  essentially  no  bioluminescent  stimulation  in  laminar  flow, 

3.  constant  peak  PMT  values  in  turbulent  flow. 

Each  of  these  observations  need  to  be  reevaluated  in  the  context  of  the  present  work. 

1.  Throughout  the  current  measurements,  an  abrupt  increase  in  bioluminescence  at  Rctrans  is 
also  often  observed  (e.g.,  figures  6a  and  7a),  particularly  when  occurring  at  high  Reynolds  num¬ 
bers.  Concerning  the  reason  for  this  increase  in  bioluminescent  stimulation  at  transition,  Rohr  et 
al.  (1990)  had  written: 
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“Whether  this  principally  results  from  the  corresponding  increase  in  shear  stress  or  the 

sudden  appearance  of  fluctuating  stresses  associated  with  the  turbulent  eddy  field 

remains  to  be  resolved.” 

Where  comparisons  can  be  made  (see  figures  7b  and  8),  the  present  PMT  measurements 
(expressed  as  counts  per  sec)  suggest  that  shear  stress,  regardless  of  whether  it  is  laminar  or  tur¬ 
bulent,  reasonably  characterizes  bioluminescent  stimuli  in  fully  developed  pipe  flow.  Therefore, 
when  plotted  as  a  function  of  Re  number,  the  increase  in  bioluminescence  observed  at  transition 
appears  to  be  at  least  partly  accounted  for  by  the  accompanying  jump  in  Vaii  (presumably  the 
increased  mixing  and  reduced  average  velocity  associated  with  turbulence  has  some  effect  also). 
For  fully  developed  pipe  flow,  Xwaii  increases  by  a  factor  of  0.005(Re)^/'^  at  transition.  The 
increase  in  shear  stress  for  Retrans  values  of  6000  and  10000  are,  respectively,  by  factors  of  3.4 
and  5.0  throughout  the  pipe. 

Another  example  of  increased  bioluminescence  occurring  at  transition  may  be  the  following. 
Dolphins  swimming  through  phosphorescent  seas  have  been  reported  to  produce  interesting  dis¬ 
plays  of  light  (Stevens,  1950;  Hill,  1950;  Thompson,  1971;  Ridgeway  and  Carter,  1993).  Pro¬ 
nounced  bioluminescence  beginning  below  the  dorsal  fin  of  a  swimming  Pacific  white-sided  dol¬ 
phin  has  been  reported  (Wood,  1973).  On  a  dolphin’s  spindle-shaped  body,  Retrans  might  be 
expected  to  occur  in  this  area  (Ridgeway  and  Carter,  1993),  and  would  result  in  both  higher 
shear  stresses  and  a  thicker  boundary  layer.  A  representative  Retrans  value  for  a  dolphin  is  6.5  x 
10^  (Newman  and  Wu,  1975).  For  comparison,  the  increase  in  shear  stress  along  a  flat  plate  at 
Retrans  =  6.5  X  10^  is  a  factor  of  6.0.  Note  even  when  x  >  20  dynes/cm^,  where  peak  levels  of 
individual  organisms  are  independent  of  x,  transition  may  still  be  associated  with  increased  bio¬ 
luminescence.  Increased  bioluminescence  could  result  from  a  larger  volume  of  flow  experienc¬ 
ing  shear  stresses  greater  than  threshold,  as  well  as  a  greater  fraction  of  the  organisms  present 
being  stimulated. 

2.  The  lack  of  significant  bioluminescence  reported  in  laminar  flow  by  Rohr  et  al.  (1990) 
deserves  comment.  If  a  value  of  Retrans  is  assumed  for  Rohr  et  al.’s  (1990)  PMT  vs  Re  data 
shown  in  their  figure  4,  then  Xwaii  values  for  fully  developed  laminar  (equation  2.6b)  and  turbu¬ 
lent  (equation  2.9b)  flow  can  be  estimated.  Based  on  Rohr  et  al.’s  (1990)  accompanying  dye 
visualization  experiments,  Retrans  ~  6000.  Included  with  the  present  San  Diego  Bay  data  in  fig¬ 
ure  6b,  are  the  2.2 1  average  PMT  levels  from  Rohr  et  al.  (1990)  as  a  function  of  the  estimated 
Xwaii  values.  When  plotted  with  logarithmic  coordinates  in  figure  6b,  a  rise  of  luminescence  with 
increasing  laminar  Xwaii  levels  is  clearly  evident  for  the  Rohr  et  al.  (1990)  data.  Although  the 
light  intensity  of  the  Rohr  et  al.  (1990)  measurements  shown  in  figure  6b  are  about  an  order  of 
magnitude  less  than  the  present  data,  the  similarity  in  general  dependence  of  bioluminescence  vs 
Xwaii  is  noteworthy.  The  dashed  line  in  figure  6b  represents  the  slope  (1.30,  R^  =  0.96)  of  the 
least-squares  power  law  regression  to  the  turbulent  Rohr  et  al.  (1990)  data.  Expressed  as  counts 
per  volume,  this  slope  would  be  0.73. 

3.  Rohr  et  al.  (1990)  have  also  reported  that  the  growing  bioluminescent  intensity,  which 
they  normally  observed  with  increasing  turbulent  pipe  flow,  did  not  appear  to  be  accompanied  by 
a  heightening  of  individual  flash  levels.  The  time  series  measurements  presented  here  offer  finer 
resolution  (0.1-s  time  averages)  than  those  (1.0-s  time  averages)  of  Rohr  et  al.  (1990).  Where 
wall  shear  stress  measurements  overlap,  Xwaii  >  40  dynes/cm^  (turbulent  flow),  the  present  results 
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show  a  similar  behavior.  Throughout  this  range  it  can  again  be  concluded  that  the  increase  in 
bioluminescence  with  Twaii  results  from  more  organisms  being  stimulated,  rather  than  an  increase 
in  individual  peak  intensity.  However,  fortwaii  <  10  dynes/cm^  (laminar  flow),  outside  the  mea¬ 
surement  range  of  Rohr  et  al.  (1990),  a  graded  response  by  individual  luminescent  organisms  is 
now  found.  The  data  in  figure  11  suggest  that  the  graded  bioluminescent  response  of  individual 
flashes  may  have  an  important  influence  on  the  form  of  the  100-s  average  PMT  data,  which  also 
exhibit  a  dramatic  change  in  form  around  Xwaii  =  10  dynes/cm^.  The  details  of  this  influence  can¬ 
not  be  inferred  from  figure  11,  however,  since  additional  information  regarding  the  fraction  of 
organisms  excited  as  a  function  of  Xwaii  is  not  known. 

4.2  RESPONSE  OF  BIOLUMINESCENT  ORGANISMS  IN  LAMINAR  FLOW 

Previously,  hydrodynamic  stimulation  of  bioluminescence  has  been  almost  exclusively  asso¬ 
ciated  with  turbulent  flow  (Anderson  et  al.,  1989;  Rohr  et  al.,  1990).  Bathyphotometers  devel¬ 
oped  for  in  situ  bioluminescent  measurements  have  always  provided  highly  turbulent  flow  as  the 
hydrodynamic  stimulus  (Seliger  et  al.,  1969;  Losee  et  al.,  1985;  Widder  et  al.,  1993).  The  present 
study  shows,  however,  that  laminar  flow  stimulates  bioluminescence.  The  question  arises  as  to 
why,  within  a  biological  context,  should  luminescent  dinoflagelates  be  excited  in  laminar  flow? 

Natural  dinoflagellate  bioluminescence  is  often  observed  in  the  presence  of  predatory  cope- 
pods  (Buskey  et  al.,  1983).  Studies  have  shown  (Esaias  &  Curl,  1972;  H.  White,  1979)  that 
marine  copepods  ingest  highly  bioluminescent  dinoflagellates  at  a  lesser  rate  than  their  nonbiolu- 
minescent  counterparts.  Furthermore,  copepods  have  also  been  found  to  exhibit  a  ‘startle’  or 
photophobic  response  to  natural  and  stimulated  dinoflagellate  bioluminescence  (Buskey  et  al., 
1983;  Buskey  &  Swift,  1985).  Since  most  luminous  flashes  are  several  orders  of  magnitude 
greater  than  the  visual  threshold  of  their  predators  (Nicol,  1978;  Lythgoe,  1979),  temporary 
blinding  (“flash  bulb  effect”)  by  the  flash  seems  likely  (Morin,  1983).  The  possibility  (Burken- 
road,  1943)  also  exists  that  luminescent  dinoglagellates  might  expose  the  copepod  feeding  on 
them  to  the  copepod ’s  predator  (“burglar  alarm  effect”). 

Taken  together,  these  observations  strongly  suggest  that  the  luminescent  response  of  dinofla¬ 
gellates  serves  as  a  defense  strategy  and  as  such,  should  be  stimulated  by  the  predatory  motions 
of  copepods.  These  motions  are  unambiguously  laminar.  A  representative  Reynolds  number  for 
copepod  motion,  even  at  a  pulse  speed  of  20  cm/s,  is  only  around  300  (Vogel,  1981),  which  is 
low  enough  to  ensure  laminar  flow  about  the  body.  Food  capture  by  copepods  is  also  thought 
(Koehl  &  Strieker,  1981)  to  be  viscously  dominated  and  laminar  in  nature.  In  retrospect,  consid¬ 
ering  both  the  natural  flow  fields  that  excite  luminescent  dinoflagellates  and  the  possible  advan¬ 
tage  that  their  bioluminescence  may  serve,  a  sensitivity  of  these  organisms  towards  shear  stresses 
induced  by  laminar  flow  is  to  be  expected. 

The  threshold  for  bioluminescent  stimulation  found  in  the  present  study  is  of  the  order  of  a 
few  dynes/cm^  (figure  9;  Xwaii  ~  1-2  dynes/cm^,  Xgvg  ~  0.5-1  dynes/cm^),  and  deserves  some 
discussion  as  to  the  possible  advantage  this  level  of  threshold  excitation  might  serve.  Clearly  for 
bioluminescence  to  be  an  effective  defense  mechanism  for  an  organism,  the  level  of  threshold 
stimulus  must  be  lower  than  the  excitation  field  caused  by  the  approach  of  its  predators.  Thresh¬ 
old  levels  must  also  be  much  higher  than  oceanic  background  levels  to  avoid  bioluminescent 
exhaustion.  The  possibility  of  continual  random  stimulation  by  ambient  motion  would  severely 
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limit  whatever  useful  purpose  bioluminescent  might  otherwise  have.  In  this  regard  Harvey 
(1929)  has  written: 

“who  can  suggest  the  use  of  light ...  to  a  protozoan,  living  at  the  surface  of  the  sea, 

blown  hither  and  thither  by  the  wind.” 

Estimates  regarding  the  magnitude  of  ocean  surface  shear  stress  that  marine  dinoflagellates 
are  exposed  to  have  been  previously  made  (Thomas  &  Gibson,  1990, 1992).  The  expected  range 
of  shear  stress  at  the  sea  surface  for  light  winds  is  about  two  orders  of  magnitude  less  than  the 
bioluminescent  threshold  levels  reported  here.  A  wind  of  3.5  m/s  (measured  10  m  above  the  sea 
surface)  would  only  produce  maximum  shear  stresses  at  the  ocean  surface  equal  to  one-tenth  that 
required  for  bioluminescent  stimulation.  Typical  shear  stress  levels  measured  throughout  the 
deep  surface  mixed  layer  are  around  0.001  dynes/cm^  (Evans,  1982),  about  three  orders  of  mag¬ 
nitude  less  than  that  necessary  for  bioluminescent  excitation. 

On  the  other  hand,  order  of  magnitude  estimates  of  the  average  shear  stress  experienced  by  a 
dinoflagellate  about  to  be  devoured  are  not  significantly  greater  than  these  same  threshold  val¬ 
ues.  The  entrainment  of  the  common  bioluminescent  dinoflagellate  Gonyaulax polyedra  (mod¬ 
eled  as  an  equivalent  35-pm  diameter  sphere;  Kamykowski,  1992),  by  water  swept  (Uq  =  3.2 
cm/s;  Koehl  &  Strickler,  1981)  towards  a  devouring  copepod  results  in  a  Reynolds  number 
around  1.1.  At  these  low  Reynolds  numbers,  Stokes  flow  is  a  reasonable  approximation  for  our 
purposes  (Bird  et  al.,  1960).  The  shear  stress  on  a  sphere  in  Stokes  flow  is  of  the  form  (White, 
1974): 


tre  =  -1.5  [pUo/R]  sin(0) . 

The  above  equation  yields,  for  our  Gonyaulax  polyedra  model,  average  and  maximum  surface 
shear  stresses  of  about  20  and  30  dynes/cm^  respectively.  A  similar  calculation  shows  that  the 
organism’s  own  propulsion  (~  0.03  cm/s,  Kamykowski  et  al.,  1992)  subjects  it  to  shear  stresses 
that  are  about  an  order  of  magnitude  less  than  threshold  values.  While  these  estimates  are 
encouraging,  it  must  be  realized  that  on  a  cellular  level  the  exact  mechanism  by  which  mechani¬ 
cal  stimulation  triggers  light  emission  remains  unclear.  Also,  the  present  experiments  do  not 
address  the  effects  of  acceleration  (conceivably  an  important  precursor  to  consumption),  and 
large  fluctuating  shear  stresses  for  stimulating  bioluminescence. 

Latz  et  al.  (1994)  provide  a  much  more  comprehensive  comparison  between  bioluminescent 
threshold  levels  and  various  oceanic  flow  fields.  Based  on  essentially  identical  threshold  values 
of  shear  stress  (determined  in  Couette  flow),  they  predict  that  surface  breaking  waves,  adult 
scombrid  fish,  and  copepod  burst  swimming  should  all  excite  bioluminescence.  All  of  these  pre¬ 
dictions  have  been  confirmed  by  observation. 

4.3  RESPONSE  OF  BIOLUMINESCENT  ORGANISMS  IN  TURBULENT  FLOW 

The  general  behavior  of  the  laminar  and  turbulent  flow  data  exhibited  in  figure  8  suggests  the 
following  results  if  two  pipe  flow  experiments  are  performed  under  identical  conditions,  except 
with  greatly  different  Retrans  values.  Under  these  conditions,  the  laminar  flow  data  of  the  high 
Retrans  experiment  should  both  bridge  the  transition  gap  of  the  low  Retrans  experiment  and  over¬ 
lap  its  turbulent  data.  In  other  words,  for  the  same  value  of  Twain  bioluminescent  levels  expressed 
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as  counts  per  second  do  not  appear  (at  least  for  range  of  overlap  studied,  20  dynes/cm^  <  Xwaii  < 
70  dynes/cm^)  particularly  sensitive  to  whether  the  fully  developed  pipe  flow  is  laminar  or  tur¬ 
bulent. 

In  a  fluid  dynamics  context,  it  is  somewhat  surprising  that  for  equal  Xwaii  values  in  laminar 
and  turbulent  pipe  flow,  average  levels  of  bioluminescence  are  nearly  the  same.  The  length  scale 
of  the  turbulent  eddies,  which  are  of  the  order  of  the  size  of  the  luminescent  organisms  (see  fig¬ 
ure  2.5,  appendix  2),  has  been  previously  believed  to  be  a  significant  factor  in  effective  stimula¬ 
tion  (Morita  &  Johnson,  1974;  New  Scientist,  1990;  Rohr  et  al.,  1990;  Widder  et  al.,  1993).  The 
lack  of  increased  sensitivity  to  small  scale  turbulence,  however,  is  consistent  with  the  luminous 
response  of  individual  organisms  as  inferred  from  0.1-s  peak  PMT  counts.  For  x^aii  >  10-20 
dynes/cm^  peak  flash  values  were  found  in  figure  11  to  have  remained  essentially  unchanged 
with  increasing  flow  stimulus  throughout  both  laminar  and  turbulent  flow.  This  behavior  sug¬ 
gests  that  whatever  additional  stimulus  accompanies  turbulent  flow,  it  is  largely  superfluous 
towards  eliciting  higher  peak  luminescent  levels  from  individual  organisms. 

There  remains  the  question;  For  the  same  shear  stress  profile  across  the  pipe,  why  doesn’t  the 
increased  mixing  associated  with  turbulent  flow  result  in  more  bioluminescence?  It  is  possible 
that  for  the  same  Xwaii,  the  increase  in  bioluminescence  due  to  turbulent  mixing  is  approximately 
compensated  by  the  associated  decrease  in  convective  velocity.  For  example  (see  figure  13),  con¬ 
sider  turbulent  pipe  flow  at  a  Reynolds  number  of  6,000.  The  necessary  laminar  mean  flow  to 
provide  the  same  shear  stress  profile  is  nearly  3.5  times  as  great  (Reiam  =  20,400).  As  shown  in 
appendix  3,  throughout  the  present  flow  regime  where  bioluminescence  is  excited,  the  contribu¬ 
tion  to  PMT  levels  (measured  as  counts  per  second)  due  to  the  convection  of  luminous  organ¬ 
isms  past  the  PMT  sensor  varies  nearly  linearly  with  Uavg- 

Conversely,  our  interpretation  of  the  response  of  bioluminescence  (expressed  as  counts  per 
second)  to  increasing  pipe  flow  will  depend  on  whether  the  flow  is  laminar  or  turbulent.  Con¬ 
sider,  for  example,  figure  8,  and  Xwaii  values  greater  than  15  dynes/cm^.  Here  bioluminescent 
levels  for  both  laminar  and  turbulent  flow  increase  approximately  linearly  with  Xwaii-  Since  in 
fully  developed,  laminar  pipe  flow  Xwaii  Ugvg  (equation  2.6b),  the  near  linear  increase  of  bio¬ 
luminescence  with  Xwaii  can  be  attributed  primarily  due  to  convection  alone.  For  fully  developed 
turbulent  pipe  flow,  however,  Xwaii  Uavg^'^^  (equation  2.11).  Consequently  Uavg  (or  Re)  in  tur¬ 
bulent  flow  increases  much  more  slowly  than  in  laminar  flow  over  the  same  range  of  Xwaii  (see 
figure  13).  Therefore,  in  turbulent  flow,  a  linear  increase  of  bioluminescence  with  x^aii  repre¬ 
sents  a  greater  increase  (consistent  with  the  enhanced  mixing)  than  what  the  associated  change  in 
convective  velocity  can  alone  provide. 

In  figure  14,  each  data  set  of  figure  8,  with  its  PMT  levels  expressed  as  counts  per  second 
(14al,  bl,  and  cl),  is  replotted  with  its  PMT  levels  expressed  as  counts  per  volume  (14a2,  b2, 
and  c2).  After  removing  the  influence  of  convection,  a  greater  increase  in  bioluminescence  with 
Xwaii  is  now  apparent  for  turbulent  flow  (compare  figures  14a  and  14b  with  figure  14c).  Note,  the 
difference  in  slope  between  what  is  shown  in  figure  14b2  (0.3;  which  was  calculated  from  mean 
velocity  measurements)  and  what  was  stated  previously  in  the  text  (0.0;  which  was  calculated 
from  equation  (2.6b)).  This  discrepancy  reflects  the  poorer  approximation  to  fully  developed 
pipe  flow  for  the  laminar,  high  Reynolds  number  data. 
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5.0  CONCLUSIONS 


“Although  in  most  ways  the  exact  manner  in  which  water  moves  is  difficult  to  perceive 
and  still  more  difficult  to  define,  as  are  the  forces  attending  such  motion,  certain  general 
features  both  of  the  forces  and  motions  stand  prominently  forth,  as  if  to  invite  or  defy 
theoretical  treatment.” 


Osborn  Reynolds  (1883) 

While  the  present  effort  provides  neither  an  adequate  identification  of  the  luminescent  organ¬ 
isms  involved  nor  a  detailed  description  of  their  interaction  with  the  flow  field,  nevertheless,  cer¬ 
tain  unmistakable  trends  in  the  data  appear  “as  if  to  invite  or  defy”  explanation.  These  trends 
include  the  repeated  observation  that  bioluminescence  is  stimulated  in  laminar  flow,  with  thresh¬ 
old  values  of  the  order  of  a  few  (0.5-3)  dynes/cm^.  With  respect  to  the  estimates  made  herein  of 
the  natural  hydrodynamic  stimulus  these  organisms  would  normally  be  subjected  to,  this  level  of 
threshold  stimulus  should  have  antipredation  value.  Another  general  feature  of  these  measure¬ 
ments  is  that  with  increasing  wall  shear  stress  a  conspicuous  break  in  the  rate  of  increase  of  aver¬ 
age  bioluminescence  is  normally  observed  between  twaii  =  10-20  dynes/cm^,  with  the  greater 
increase  occurring  initially.  This  break,  which  occurred  in  laminar  flow,  is  believed  to  be  at  least 
partly  related  to  the  behavior  of  the  peak  PMT  levels  associated  with  individual  bioluminescent 
flashes.  Peak  PMT  levels  initially  exhibited  a  graded  response,  increasing  with  Xwaii  from  thresh¬ 
old  to  Twaii  values  of  about  10  dynes/cm^,  and  then  remaining  nearly  constant  thereafter.  Finally, 
throughout  the  figures  for  Xwaii  >  20  dynes/cm^,  an  approximate  linear  dependence  between 
average  PMT  levels  (expressed  as  counts  per  sec)  and  wall  shear  stress  is  usually  observed, 
regardless  of  whether  the  flow  is  laminar  or  turbulent.  The  interpretation  of  the  stimulatory 
nature  of  this  relationship,  as  previously  discussed,  depends  on  the  relationship  between  Xwaii  and 
Uavg  which  is  different  in  laminar  and  turbulent  flows. 

The  few  bioassays  complementing  the  present  data  set,  as  well  as  the  historic  record,  suggest 
that  the  luminous  organisms  responsible  for  the  stimulated  light  belong  only  to  a  few  species  of 
dinoflagellates.  In  coastal  environments,  the  proportion  of  bioluminescent  species  is  only  1  to 
2%  (Morin,  1983)  of  that  found  in  oceanic  waters.  However,  higher  productivity  in  coastal 
waters  may  result  in  a  greater  amount  of  bioluminescence  potential  per  unit  volume  of  near  sur¬ 
face  water  (Buskey  et  al.,  1987).  In  reference  to  figures  6b,  7b,  8,  and  11  it  is  compelling  to  sug¬ 
gest  that  there  may  be  certain  regimes  of  flow  agitation  where  coastal  dinoflagellates  share  simi¬ 
lar  bioluminescent  fingerprints  (i.e.,  similar  peak  and  average  PMT  dependance  on  x^aii)-  It  is 
hoped  that  this  work  will  lead  to  better  predictive  capabilities  towards  assessing  bioluminescent 
stimulation  in  a  variety  of  oceanic  flow  fields.  Further  work,  particularly  with  offshore  organ¬ 
isms  and  cultures  of  luminous  dinoflagellates  (where  species,  concentration,  and  prestimulus  his¬ 
tory  can  be  better  controlled)  is  necessary  to  corroborate  this  initial  study. 
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Figure  1 .  Photograph  of  bioliiminescence.  The  marine  bioluminescent  dinoflagellate,  Gonyaulax  polyedra  (culture  of  ~10®  cells  per  liter),  are  stimulated  as  they 
poured  over  a  17.5  cm  model  of  a  dolphin. 


a.  The  outline  of  a  rope. 


b.  A  0.25-m  diameter  sphere. 


c.  A  5-gallon  box.  d.  The  presence  of  divers. 

Figure  2.  Computer  enhanced  photographs  of  in  situ,  flow  induced  bioluminescence  obtained  via  a  low-light  video 
camera.  Objects  can  be  identified  3-5  meters  beneath  the  sea  surface  as  a  result  of  stimulated  bioluminescence.  At 
these  depths  it  is  likely  that  one  or  more  species  of  dinoflagellates  are  responsible  for  the  emitted  light. 


Figure  3.  Diagrammatic  illustration  of  the  primary  pipe  flow  apparatus 


Reynolds  Number 

Figure  4.  Representative  seawater  friction  factor  vj  Reynolds  number  data.  Most  of  the  measurements  show 
good  agreement  with  theoretical  (laminar  flow  -  solid  line)  and  empirically  derived  (turbulent  flow  -  dashed 
line)  values.  Laminar  flow  measurements  taken  in  the  0.635-cm  i.d.  pipe  ( •  ,  ♦  ,  A )  for  Re  >  2000,  howev¬ 
er,  begin  to  diverge  from  theoretical  values  because  the  pipe’s  dimensionless  entrance  length  (L(./D)  is  not 
long  enough  to  support  fully  developed  flow.  Measurements  obtained  in  the  0.305-cm  i.d.  pipe  ( ▼),  with  its 
much  longer  entrance  length,  exhibit  better  agreement. 


Reynolds  Number 


Figure  5.  Downstream  flash  termination  location  of  bioluminescent  organisms  stimulated  at  the  pipe  inlet  as  a 
function  of  Reynolds  number.  A  response  time  of  100  ms  is  assumed,  followed  by  flash  durations  of  100  ms 
(short  dash)  and  200  ms  (long  dash).  PMT  positions  are  indicated  by  solid  lines.  Pipe  i.d.  =  0.635  cm;  69°F 
(v  =  0.0104  cm^/sec) 
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Figure  6a.  Representative  San  Diego  Bay  PMT  levels  vi  Reynolds  number.  The  PMT  is  located  67  cm  from 
the  inlet  of  the  0.635-cm  i.d.  pipe.  Each  data  point  represents  an  average  of  2.2-1  of  seawater.  Open  symbols 
(Run  1;  □ ,  A)  refer  to  measurements  taken  from  the  bottom  half  of  the  75-1  tank,  where  the  flow  sequence  is 
monotonically  decreasing.  Solid  symbols  (Run  2;  ■  ,  A )  refer  to  measurements  taken  from  the  remaining 
half  of  the  tank,  where  the  flow  sequence  is  monotonically  increasing.  Transition  (dotted  line)  from  laminar  to 
turbulent  flow  occurs  at  Retr^ns  =  6000,  where  a  change  in  the  PMT  vj  Re  dependence  is  evident. 
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Figure  6b.  PMT  levels  wall  shear  stress.  San  Diego  Bay  data,  collected  67  cm  from  the  0.635-cm  i.d. 
pipe  inlet,  are  represented  by  symbols  and  composed  both  of  2.2-1  (  □ ,  A;  same  PMT  data  as  in  figure 
6a)  and  100-s  (♦,©,■)  averages.  Sea  of  Cortez  data,  collected  65  cm  from  the  0.635-cm  i.d.  pipe 
inlet,  are  represented  by  crosses  and  composed  of  2.2-1  averages.  Transition  (dotted  line)  from  laminar 
to  turbulent  flow,  for  all  the  experiments,  occurs  around  a  Reynolds  number  of  6000. 
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Figure  7a.  Representative  San  Clemente  Island  PMT  levels  vi  Reynolds  number.  Each  data  point  represents  a 
100-s  average  of  two  PMTs,  located  22  and  67  cm  from  the  inlet  of  the  0.635-cm  i.d.  pipe.  Transition  from 
laminar  to  turbulent  flow  varies  due  to  ship  motion;  for  the  open  symbol  data  ( □  ,  A,  V)  Rctran,  =  3500,  for 
the  cross  symbol  data  Rej^^nj  =  5400,  and  for  the  solid  symbol  data  ( A.  ▼)  Re,ran.s  ~  7600.  Note  that  the  PMT 
levels  for  the  high  laminar  flow  Reynolds  number  data  ( A ,  ▼ )  are  conspicuously  low. 
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Figure  7b.  PMT  levels  (same  data  as  in  figure  7a)  vs  wall  shear  stress.  The  solid  line  is  taken  from  figure  6b 
and  is  representative  of  San  Diego  Bay  data.  The  disparity  in  PMT  levels  between  laminar  and  turbulent 
flow,  evident  in  figure  7a  when  plotted  as  a  function  of  Re,  is  no  longer  apparent. 
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Figure  7c.  PMT  levels  vj  wall  shear  stress  (same  data  as  in  figure  7b)  with  laminar  (  □),  turbulent  (  ■),  and  transitional  (  flow  indicated. 


Figure  8.  PMT  levels  vs  wall  shear  stress  (x^an)  for  the  highest  laminar  pipe  flow  stimulus.  Solid  symbols 
designate  laminar  flow  measurements  taken  in  San  Diego  Bay  with  the  0.305-cm  i.d.  pipe  (  ▼,  500-s 
averages,  PMT  is  106  cm  downstream  of  inlet)  and  the  0.635-cm  i.d.  pipe  (  ■  ,  2.2-1  averages,  PMT  is 
67  cm  downstream  of  inlet),  and  off  San  Clemente  Island  with  the  0.635-cm  i.d.  pipe  (  • ,  100-s  aver¬ 
ages,  PMT  67  cm  is  downstream  of  inlet).  Open  symbols  designate  turbulent  flow  measurements  taken  in 
San  Diego  Bay  with  the  0.635-cm  i.d.  pipe  (  □ ,  2.2-1  averages,  PMT  is  67  cm  downstream  of  inlet). 

Note  that  for  >  15  dynes/cm^,  the  PMT  vs  dependence  is  similar  (approximately  linear)  regard¬ 
less  of  whether  the  flow  is  laminar  or  turbulent. 
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Figure  9.  PMT  levels  (100-s  averages)  vj  wall  shear  stress  (t„,|||)  for  the  lowest  laminar  pipe  flow  stimulus. 
San  Diego  Bay  water  is  used  in  the  0.635-cm  i.d.  pipe-flow  apparatus  with  the  PMT  positioned  67  cm  from 
the  inlet.  The  solid  line  represents  a  least  squares  power  law  regression  to  the  data  between  0.7  and  14.0 
dynes/cm2  and  has  a  slope  of  2.57.  Threshold  values  are  seen  to  be  around  1  dyne/cm^. 
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Figure  10.  Sample  laminar  (a)  and  turbulent  (b)  PMT  time  series  data  (0.1-s  bins).  The  PMT  is  located  67  cm 
from  the  inlet  of  the  0.635-cm  i.d.  pipe.  Each  vertical  spike  is  believed  to  result  from  an  individual  biolumi- 
nescent  flash.  Corresponding  wall  shear  stress  values  are  15  dynes/cm^  (a)  and  101  dynes/cm^  (b).  It  is  note¬ 
worthy  that  although  the  average  bioluminescent  level  for  (b)  is  nearly  five  times  larger  than  (a),  peak  values 
remain  about  the  same. 
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Figure  11.  Average  (100  s)  peak  (0, 1  s)  PMT  levels  as  a  function  of  wall  shear  stress  (x^an).  The 
PMT  is  located  67  cm  from  the  inlet  of  the  0,635-cm  i.d.  pipe.  Solid  symbols  (  ^ )  designate  average 
PMT  levels  recorded  in  San  Diego  Bay,  open  symbols  (0.0)  designate  peak  PMT  levels  recorded 
in  San  Diego  Bay,  and  cross  symbols  designate  peak  levels  recorded  off  San  Clemente  Island.  Note 
that  the  leveling  off  of  peak  PMT  levels  beginning  at  =  10-20  dynes/cm^,  coincides  with  the 
reduced  slope  of  the  average  PMT  vs  x„^]|  data. 
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Figure  12.  Comparison  of  peak  (0. 1  s  average)  PMT  V5  Reynolds  number  dependence  for  top  (6  cm  from 
inlet)  and  bottom  (67  cm  from  inlet)  PMT  locations  on  the  0.635-cm  i.d.  pipe. 
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Figure  13.  Reynolds  number  v^  wall  shear  stress  for  fully  developed  laminar  (solid  line)  and  turbulent  (dashed  line)  flow  in  a  0.635-cm  i.d. 
pipe,  with  seawater  at  20.6°C  (U^vg  =  0.0164  x  Re). 
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Figure  14a.  Laminar  flow  data  are  from  San  Diego  Bay. 


Figure  14.  Comparison  of  PMT  levels  expressed  as  counts  per  second  (same  data  as  in  figure  8)  and  counts  per  vol 
ume  as  a  function  of  wall  shear  stress.  Laminar  flow  data  are  from  San  Diego  Bay  (14a)  and  off  San  Clemente 
Island  (14b).  Turbulent  flow  data  are  from  San  Diego  Bay  (14c). 
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Figure  14b.  Laminar  flow  data  off  San  Clemente  Island. 
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Figure  14c.  Turbulent  flow  data  are  from  San  Diego  Bay. 
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APPENDIX  1:  VIDEO  CALIBRATION  AND  LIGHT  ESTIMATES  FOR 
PROJECT  BONS  (Bioluminescence  for  Ocean  Night  Surveillance) 


1.1  INTRODUCTION 

NRaD,  Code  524,  and  the  Naval  Oceanographic  Office  (NAVOCEANO),  Code  N3212,  have 
been  involved  in  bioluminescence  measurements  in  many  areas  of  the  Northern  Hemisphere  (ref¬ 
erences  3-17).  This  work  shows  that  night  bioluminescence  is  omnipresent  in  the  marine  envi¬ 
ronment,  and  most  intense  in  near  shore  regions  and  in  the  upper  mixed  layer  of  the  ocean  (upper 
20-50  meters)  (references  3,  4,  6,  7,  9, 12-17).  Luminous  plankton  are  known  to  be  easily  stimu¬ 
lated  (references  1,  6,  8, 10)  and  stimulated  bioluminescence  from  swell  and/or  current-induced 
flow  has  been  observed  around  submerged  objects  (references  6,  18-21). ^  An  ISIT  video  camera 
was  used  successfully  to  image  bioluminescence  generated  from  submerged  objects  at  night  dur¬ 
ing  tests  at  sea  (reference  2).^  Four  at-sea  operations  have  been  completed  and  video  pictures  of 
stimulated  bioluminescence  around  objects  were  obtained.  This  appendix  presents  the  procedure 
leading  to  the  calibration  curves  which  allow  in  situ  bioluminescence  levels  to  be  estimated.  As 
an  example,  a  bioluminescence  image  of  a  sphere  taken  from  the  cruise  of  February /March  1989 
in  the  Gulf  of  California  is  analyzed. 

Methods 

A  Cohu  Inc.,  ISIT  camera,  model  #5162-2000/0000  and  serial  #144104,  and  a  Hamamatsu 
video  analyzer,  model  #DVS-3000  and  serial  #88020041,  were  obtained  as  part  of  NRaD  Code 
524  capital  equipment.  The  camera  has  both  automatic  and  manual  controls  for  both  the  video 
gain  and  intensifier  high  voltage.  The  video  gain  and  high  voltage  settings  were  digitally  read  so 
camera  settings  can  be  reproduced  to  allow  for  calibration.  The  video  and  intensifier  gains  were 
manual  and  set  at:  video  +0.02  V  and  HV  +0.04  V.  The  DVS-3000  video  analyzer  brightness 
was  set  at  0.5  (center  of  range)  and  the  contrast/gain  was  set  at  1.0  (center  of  range).  Two  Fuji- 
non  lenses  were  used  on  the  camera:  50  mm  f/0.7  and  25  mm  f/0.85,  aperture  areas  40  and  6.8 
cm^  respectively. 

The  DVS-3000  was  used  to  put  date/time  information  on  each  frame  during  data  acquisition. 
During  video  playback,  the  analyzer  captures  (digitizes)  individual  frames  and  allows  pixel  data 
from  areas  within  a  video  frame  to  be  integrated.  The  integration  allows  a  quantified  video  signal 
to  be  obtained  which  can  be  calibrated  to  an  absolute  light  level. 

Quantifiable  light  sources  were  prepared  from  solutions  of  luminescent  marine  bacteria.  One 
ml  samples  of  the  solution  were  measured  in  a  Quantalum  luminescence  photometer  calibration 
unit  (model  #2000  and  serial  #882).  The  Quantalum  unit  itself  was  calibrated  by  using  secondary 
light  standards  (Ni^^  in  conjunction  with  phosphors).  The  one  ml  bacteria  solutions  measured  in 
the  Quantalum  calibration  unit  provide  a  known  source  strength  (photons/sec).  The  one  ml  bac¬ 
teria  solutions  were  placed  at  a  distance  of  5.5  m  from  the  ISIT  camera  lens  and  the  image  was 
recorded  using  a  VCR. 

'  For  induced  laminar  flows  of  as  little  as  1/2  kt  (26  cm/s),  wall  shear  stress  on  a  20-cm-long  flat  plate  will 
everywhere  exceed  bioluminescence  threshold  levels  (—1  dyne/cm^,  as  determined  in  the  work  herein). 

2  NOSC  Technical  Notes  (TNs)  are  working  documents  and  do  not  represent  an  official  policy  statement  of 
NCCOSC  RDT&E  Div.  (formerly  NOSC).  For  further  information,  contact  the  author. 
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To  quantify  the  video  signals  and  the  corresponding  variance,  800  frames  of  the  imaged 
bacteria  solution  were  averaged  100  frames  at  a  time.  The  image  spot  (20  x  20  pixels)  of  each 
100  frame  average  was  integrated.  A  similar  20  x  20  pixel  area  was  integrated  without  a  source 
to  determine  the  background  video  signal.  The  source  video  signals,  minus  background  of  these 
eight  integrations,  were  then  averaged  to  obtain  the  final  video  signal  for  each  bacteria  light 
source.  The  bacteria  source  was  measured  immediately  before  and  after  the  image  was  recorded, 
and  these  Quantalum  readings  were  averaged  to  obtain  the  source  strength. 

The  one  ml  bacteria  solution  source  strength  was  divided  by  4jtr^  and  multiplied  by  the  lens 
area  to  obtain  the  solid  angle,  and  thus  the  light  entering  the  lens  (photons/sec).  The  number  of 
pixels  illuminated  in  the  image  is  four  times  greater  for  the  50  mm  lens  relative  to  that  of  the 
25  mm  lens.  This  difference  in  pixels  is  due  to  the  x2  magnification  of  the  50  mm  lens  relative 
to  the  25  mm  lens.  The  number  of  pixels  for  the  images  was  estimated  to  be  100  (10  x  10)  for  the 
25  mm  lens  and  400  (20  x  20)  for  the  50  mm  lens.  To  express  the  calibration  in  fundamental 
units,  the  light  entering  the  lens  and  the  video  signal  were  divided  by  the  estimated  image  pixels 
to  give:  photons/sec/pixel  vs  video  signal  per  pixel.  A  range  of  light  intensities  was  obtained  by 
using  various  dilutions  of  the  bacteria  solutions  and  the  resulting  calibration  curves  (Video  ISIT 
calibration  50  mm  lens,  figure  1.1a  and  Video  ISIT  calibration  25  mm  lens,  figure  1.1b)  were 
obtained.  A  second  order  polynomial  curve  was  fit  to  each  calibration  curve  as  shown  in  figures 
1.1a  and  1.1b. 

The  camera  video  calibration  provided  an  estimate  of  the  bioluminescence  signal  generated 
by  a  24-cm  diameter  black  float.  Video  pictures  were  obtained  on  the  nights  of  2/26/89  and 
2/27/89  in  the  southern  Gulf  of  California  about  40  miles  off  La  Paz,  Mexico.  There  was  back¬ 
ground  light  from  the  upper  levels  of  the  ship  (USNS  De  Steiguer),  the  sky  was  clear  with  no 
moon,  and  there  was  little  or  no  wind.  The  surface  water  was  smooth  with  a  slight  swell  and  the 
depth  was  on  the  order  of  1000  m.  The  bioluminescence  measured  (about  1E8  photons/sec/ml) 
with  the  NRaD  Code  524  photometer  system  was  typical  for  the  area.  The  float  was  a  Benthos 
glass  sphere  wrapped  in  black  shade  cloth  and  susprended  from  a  black  electrical  cable  (1.4  cm  in 
diameter)  with  a  weight  suspended  below  it.  The  float  was  lowered  from  the  ship’s  side  winch. 
The  video  camera  was  mounted  on  a  tripod  which  was  on  the  second  deck  at  a  height  of  5.3  m 
above  the  sea  surface  and  pointed  down  at  about  30  degrees  from  the  vertical.  The  float  was  sus¬ 
pended  at  a  depth  of  6.1  m  and  at  a  horizontal  distance  of  about  3.6  m  from  the  camera. 

A  video  frame  containing  the  float  image  was  captured  (figure  2c)  and  a  50  x  50  pixel  square 
centered  around  the  image  was  integrated.  A  similar  50  x  50  pixel  square  in  the  background  was 
integrated  to  obtain  an  equivalent  background  signal.  The  background  signal  was  subtracted 
from  that  of  the  float  image  to  give  a  net  video  signal  of  about  9E4.  The  average  video  signal  per 
pixel  is  then  36,  which  is  equivalent  to  175  photons/sec/pixel  from  the  50  mm  lens  calibration 
curve.  The  light  flux  at  the  lens  is  then  175  ph/s/pixel  *  2500  pixels/40  cm^  aperture  which 
equals:  1.09E4photons/sec/cm^. 

The  next  step  is  to  determine  what  fraction  of  the  bioluminescence  from  the  submerged 
object  gets  to  the  camera  lens.  This  was  done  by  writing  a  Quick  Basic  program  (the  program 
VIDEOATT.BAS  listing  is  included).  The  program  calculates  the  primary  ray  light  path  from  the 
submerged  object  to  the  camera  lens  and  then  calculates  the  light  attenuation  along  that  path. 
Having  the  light  path,  the  attenuation  fraction  is  determined  by  the  following  calculations: 
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1.  The  fraction  of  light  transmitted  through  the  water/air  interface.  For  this  case:  0.95 

2.  The  inverse  square  distance  from  the  underwater  virtual  image  to  the  lens.  For  this  case: 
7.38E-8/cm2 

3.  The  attenuation  due  to  the  water  path  (the  attenuation  length  is  assumed  to  be  10  m  which 
is  typical  for  the  area).  For  this  case:  0.53 

The  bioluminescence  estimate  is  then  simply  the  light  flux  at  the  camera  lens  divided  by  the 
attenuation  fraction,  which  in  this  case  is  [(1.09E4  photons/sec)/cm^]/[3.74e-8/cm^],  which 
gives  2.9ell  photons/sec.  The  bioluminescence  can  be  estimated  from  other  video  images  in  the 
same  way. 
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Figure  1.1a.  Video  ISIT  camera  calibration  curve  for  50  mm  lens. 


ph/s/pix  =  44.55  +  2.88*vs/pix  +  0.01  1  28*(vs/ plx)2 


' Program ;  VIDEOATT . BAS 
CLEAR  :  CLS 

prog$  =  "Program  to  calculate  source  strengths  underwater 
10/12/94" 

PRINT  prog$ 
restart: 

H  =  17.34  'height  above  sea  surface  in  feet:  89  Baja 
XO  =  12  'horizontal  distance  lens  to  object  in  feet:  89 

Baja 

INPUT  "Attenuation  length  in  feet  for  light  (480  nm)  in 
seawater  <33ft  (10m)>";  attlen 

IF  attlen  =  0  THEN  attlen  =  33 
'default  attenuation  length 

INPUT  "What  is  the  depth  of  the  object  in  feet";  D2 


'determines  the  actual  horizontal  distance  (X)  the  primary 
ray  leaves  water  surface 

'this  distance  (X)  determines  the  angle  of  the  primary  ray 
due  to  refraction  above  the  water  surface 


HS  =  H  2 
DS  =  D2  "  2 
wis  =  (4  /  3)  "  2 
FOR  i%  =  1  TO  500 
X  =  .01  *  i%  +  XO  -  8 
actual  horizontal  distance 
XOS  =  (XO  -  X)  "2 
at  surface  from  object  squared 
XS  =  X  "  2 
at  surface  squared 

E  =  DS  *  XS  +  XS  *  XOS  -  w 
in:  index^2  —  ( sin ( refracted) / 


'height  squared 
' depth  squared 
'water  index  squared 
'iteration  loop 

'start  iteration  8  ft  less  than 

'horizontal  distance  primary  ray 

'horizontal  distance  primary  ray 

s  *  XOS  *  (HS  +  XS)  'error 

in ( incident )) ^2  =  0 


E  =  ABS(E) 

'absolute  value  of  error 
IF  i%  =  1  THEN  EM  =  E 


'start  at  initial 


error 

IF  E  <  EM  THEN  im%  =  i% :  EM  =  E  'find  minimum  error 

NEXT  i%  'end  of  iteration  loop 


PRINT  "i%  at  minimum  =  im%;  "  error=  EM  'print 
loop  count  and  minimum  error 

X  =  .01  *  im%  +  XO  —  8  'calculate  X  from 


loop  count  at  minimum  error 

PRINT  "X=  X;  "ft",  X  /  3.28;  "m"  'print  X  in  ft 

and  metres 

phi  =  ATN(X  /  H)  'refracted  angle 

in  radians  from  X  and  H 

PRINT  "phi=  ";  57.296  *  phi  'print  refracted 

angle  in  degrees 


GOSUB  intensity 
fractional  intensity 
INPUT  q$ 
restarting 


'gosub  to  calculate 
' input  dummy  before 
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'go  back  to  start  to 


GOTO  restart 
calculate  another 
intensity: 

'fraction  transmitted  through  water  surface: 

waterindex  =4/3  'index  for  water 

=  4/3 

SINincident  =  SIN (phi)  /  waterindex  'index 

=  sin (refracted) /sin (incident) 

cosine  =  SQR(1  -  SINincident  "  2)  'cos  = 

sqr  ( 1— sin''2  ) 

COSref  =  SQR(1  -  SIN(phi)  2)  'cos  = 

sqr  ( 1— sin''2  ) 

RN  =  COSref  /  COSinc  'ratio  of 

cos ( refracted ) /cos ( incident ) 

TS  =  (4  /  SQR(3))  /  (RN  +  waterindex)  *  (RN) 

'transmitted  wave  amplitude  polarization  perpendicular  to  plane 
of  rays 

TP  =  (4  /  SQR(3))  /  (RN  *  waterindex  +  1)  *  (RN) 

'transmitted  wave  amplitude  polarization  in  plane  of  rays 
fractrans  =(TS"2+TP"'2)/2 
'intensity  transmitted 

PRINT  "fraction  intensity  transmitted  through  water  surface= 
fractrans 

'inverse  square  plus  water  attenuation: 

PRINT 

pathuw  =  SQR((X0  -X)"2+D2"2)  'actual  path 

length  underwater 

auw  =  EXP (—pathuw  /  attlen)  'attenuation 

underwater 

PRINT  "underwater  attenuation  factor=  auw 

invSqrpath  =  1  /  ( ( (XO  -  X)  /  SIN(phi))  +  (SQR(X  "  2  +  HS))) 

2  '  inverse  square  ratio  from  virtual  source  to  camera 

FractionAtLensFt  =  auw  *  invSqrpath  *  fractrans  *  1  /  (4 

*  3.14)  'per  sq  ft 
PRINT 

FractionAtLens  =  FractionAtLensFt  /  (144  *  2.54  2) 

'per  sq  cm 

PRINT  "Fractional  Light  Flux  from  object  at  lens= 
FractionAtLens;  "  /cm''2" 

PRINT 

INPUT  "Input  photons/sec/cm''2  into  camera  lens  from 
video  signal  calibration  of  image";  VSphotons 

PRINT  "Source  strength  from  object  =  ";  VSphotons 

/  FractionAtLens;  "  ph/s" 

RETURN 

END 
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APPENDIX  2:  PROPERTIES  OF  FULLY  DEVELOPED  PIPE  FLOW 


2.1  FULLY  DEVELOPED  FLOW 

Pipe  flow  is  considered  fully  developed  when  its  velocity  profile  is  no  longer  evolving  down¬ 
stream.  With  no  mean  acceleration,  the  pressure  drop,  dP/dx,  across  the  face  of  any  concentric, 
cylindrical  control  volume  of  radius  r  must  be  balanced  by  the  shear  stress  x  acting  on  its  perime¬ 
ter.  This  relationship  can  be  expressed  as 

T(r)  =  (dP/dx)  *  (r/2),  (2.1a) 

which  at  the  pipe  wall,  R,  becomes 

■Cwall  =  (dP/dx)  *  (R/2).  (2.1b) 

Eliminating  dP/dx,  equations  (2.1a)  and  (2.1b)  reduce  to 

T(r)  =  t^aii  (r/R).  (2.1c) 

This  linear  stress  distribution  expresses  a  general  law,  regardless  of  whether  the  flow  is  laminar 
or  turbulent. 

Fully  developed  pipe  flows  are  approached  asymptotically  in  the  downstream  direction,  and 
experiments  have  determined  appropriate  entrance  lengths  (Lg).  The  accepted  relations  for  lami¬ 
nar  and  turbulent  pipe  flows  are  the  following  (F.  White,  1979): 

(Le/D)iaminar  0.06  Re,  (2.2a) 

and 


(Le/D)turbulent  ~  4.4  (Re)^^^,  (2.2b) 

where  D  is  the  pipe  diameter  and  Re  is  the  Reynolds  number.  The  Reynolds  number  is  indicative 
of  the  ratio  of  the  force  acting  on  a  small  parcel  of  fluid,  to  the  force  developed  on  the  parcel  by 
the  viscous  stresses,  and  has  become  the  basic  dynamic  characteristic,  qualifying  the  state  of 
fluid  motion  (Childress,  1981).  Here  Re  is  defined  as 

Re  =  pUavgD/|x,  (2.3) 

where  p  is  the  fluid  density,  Ugvg  is  the  average  flow  speed,  and  p  is  the  kinematic  viscosity. 

2.2  LAMINAR  FLOW 

The  sole  source  of  shear  stress  in  laminar  flow  is  due  to  viscosity.  The  viscous  stresses  devel¬ 
oped  in  laminar  flows  of  most  common  fluids  (such  as  air,  water,  and  sea  water)  are  found  to  be 
equal  to  the  product  of  the  kinematic  viscosity  and  the  local  velocity  gradient.  This  relationship 
was  first  proposed  by  Isaac  Newton  in  1687  (Cajori,  1946): 

“The  resistance  arising  from  want  of  lubricity  in  the  parts  of  the  fluid,  is,  other  things 
being  equal,  proportional  to  the  velocity  with  which  the  parts  of  the  fluid  are  separated 
from  one  another.” 
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Fluids  adhering  to  this  linear  relationship  between  shear  stress  and  velocity  gradient^  are  referred 
to  as  Newtonian.  In  fully  developed  laminar  pipe  flow  of  a  Newtonian  fluid  the  following  rela¬ 
tionship  holds: 

t^laminar  ~  T:vvall(r/R)  ~  p(du/dr),  (2-4) 

where  u(r)  is  the  velocity  along  the  pipe  axis,  measured  at  radius  r  from  the  pipe  centerline.  (As  a 
matter  of  notation,  u(r),  U(r),  and  u'(r)  refer,  respectively,  to  the  instantaneous,  average,  and  fluc¬ 
tuating  velocity  at  r.  The  average  velocity  through  the  pipe  is  designated  Uavg-)  For  fully  devel¬ 
oped  laminar  flow,  equation  (2.4)  predicts  (Schlichting,  1979)  both  the  correct  parabolic  velocity 
profile  and  dimensionless  relationship  between  pressure  drop  (Darcy  friction  factor,  K)  and  aver¬ 
age  flow  rate  (Reynolds  number.  Re).  The  Darcy  friction  factor  is  defined^  as 

X  =  (dp/dx)D/(l/2pUavg2).  (2.5) 

The  theoretical  relation  between  X  and  Re  for  fully  developed  laminar  pipe  flow  is 

X  =  64/Re  (2.6a) 

which  can  also  be  expressed,  via  equations  (2.1b),  (2.5),  and  (2.6a),  as 

twall  lam  ~  SpUavg^/Rc  Uavg-  (2.6b) 


2.3  TURBULENT  FLOW 

In  turbulent  flow,  another  mechanism  for  momentum  exchange  arises.  Turbulent  shear  stress^ 
results  from  the  correlation  of  two  components  of  velocity  fluctuation  at  the  same  point.  For 
fully  developed  turbulent  pipe  flow  the  total  shear  stress  becomes 

ttotai  =  VaiKr/R)  =  p(du/dr)  +  p  u'v'  (2.7) 

where  u  v  represents  the  time-averaged  product  of  the  fluctuating  velocities  in  the  longitudinal 
and  radial  directions.  While  individual  velocity  fluctuations  are  only  a  few  percent  of  the  average 
mean  velocity"^,  turbulent  shear  stresses,  except  within  the  viscous  sublayer  at  the  wall,  essen¬ 
tially  constitute  the  total  shear  stress  throughout  the  pipe  (Laufer,  1954).  This  is  illustrated  in  fig¬ 
ure  2.1  where  the  theoretical,  total  shear  stress  (solid  line),  identical  for  both  laminar  and  turbu¬ 
lent  flow,  is  compared  with  the  turbulent  shear  stress  measurements  of  Laufer  (1954).  Both  the 
total  and  turbulent  shear  stress  shown  in  figure  2.1  are  normalized  by  Xwaii-  The  Reynolds  num¬ 
bers  for  Laufer’s  (1954)  pipe  flow  experiments  were  50,000  and  500,000.  Note  that  on  the  pipe 
wall  the  velocity  is  zero  (no  slip  condition).  Therefore,  as  the  wall  is  approached,  equation  (2.7) 
will  reduce  to  equation  (2.4). 

^  In  the  limit  of  infinitesimal  changes  the  shear  strain  rate  is  equal  to  the  velocity  gradient  (White,  1979). 

2  Unfortunately,  a  dimensionless  shear  stress  is  sometimes  used  instead  of  X.  This  dimensionless  shear  stress,  f,  is 
referred  to  as  the  Fanning  friction  factor  and  is  equivalent,  via  equation  lb,  to  1/4  X. 

^  The  correlation  term  u  v  is  commonly  called  a  turbulent  stress  because  it  has  the  same  dimensions  as  the 
laminar  stress  and  has  the  same  mathematical  effect.  Actually  u  v  is  a  convective  acceleration  term  which  is  why 
the  density  appears  (White,  1979). 

The  peak  longitudinal,  radial,  and  azimuthal  turbulent  velocity  fluctuations  occur  near  the  pipe  wall  (r/D  ~ 

0.001),  are  respectively  about  8%,  4%,  and  6%  of  the  centerline  velocity,  and  decrease  monotonically  to  zero  at  the 
wall  and  to  about  3%  at  the  pipe  centerline  (Laufer,  1954). 
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The  nature  of  turbulent  shear  stress  greatly  increases  mixing.  Whereas  in  laminar  flow,  fila¬ 
ments  of  dye  remain  intact  as  they  pass  down  the  pipe,  at  the  onset  of  turbulence  they  appear  to 
“explode”,  mixing  rapidly  across  it.  Figures  2.2a  through  c  are  photographs  of  dye  taken  30  pipe 
diameters  (D  =  0.635-cm)  from  the  pipe  inlet  where  the  dye  was  introduced  for  respective 
Reynolds  numbers  of  13,500  (beginning  of  turbulent  flow),  10,500  (beginning  of  transitional 
flow),  and  3200  (laminar  flow).  At  higher  Reynolds  numbers,  and  at  the  same  position  along  the 
tube,  the  dye  streak  could  no  longer  be  distinguished.  At  the  pipe  entrance  the  filament  was 
always  intact^,  regardless  of  the  Reynolds  number.  Toward  the  end  of  the  pipe  the  dye  streak  was 
always  distinguishable  for  laminar  flows,  intermittently  intact  for  transitional  flows,  and  com¬ 
pletely  mixed  for  turbulent  flows. 

The  photographs  in  figures  2.2a  through  c  were  taken  18  hours  after  the  water  had  been 
transferred  to  the  aquarium.  Normally  experiments  began  45-60  minutes  after  the  aquarium  was 
filled  and  transition  was  usually  observed  to  occur  at  the  significantly  lower  Reynolds  number  of 
6000.  Lower  Reynolds  number  transition  resulting  from  the  residual  motion  associated  with  the 
filling  of  the  reservoir  is  a  common  occurrence  (Reynolds,  1883).  Reynolds  (1883)  had 
achieved,  after  painstaking  effort,  transitional  Reynolds  around  12,000  (Prandtl,  1934;  Rouse, 
1938).  It  is  interesting  to  note  that  when  Reynolds’  experiment  was  repeated  using  the  original 
apparatus  a  century  later,  such  relatively  high  transition  values  could  not  be  duplicated  due  to  the 
increased  background  vibrations  introduced  by  modern  traffic  (Van  Dyke,  1982)^. 

Another  fundamental  difference  between  turbulent  and  laminar  shear  stress  is  that  the  former 
is  not  determined  by  the  local  mean  velocity  gradient,  but  now  is  primarily  dependent  on  the  pre¬ 
vious  history  of  the  flow  which  carries  the  turbulent  eddies.  Therefore  mean  and  turbulent  veloc¬ 
ity  profiles,  and  the  functional  relationship  between  wall  shear  stress  and  turbulent  flow  rate, 
must  be  empirically  determined.  It  has  been  found  (Schlicting,  1979)  for  the  range  of  Reynolds 
numbers  important  to  this  study.  Re  <  20,000,  that  a  suitable  representation  for  the  mean  velocity 
profile  (except  very  near  the  wall)  for  fully  developed  turbulent  pipe  flow  is  the  following: 

U(r)/Un,ax  =  (r/R)i/^  (2.9) 


where  Umax  is  the  centerline  velocity. 

Figure  2.3a  compares  the  theoretical  laminar  and  experimentally  derived  turbulent  pipe 
velocity  profiles,  normalized  by  the  centerline  velocity,  with  the  measurements  of  Patel  and 
Head  (1969)  at  Reynolds  numbers  of  2015  and  4060.  Figure  2.3b  illustrates  the  laminar  and  tur¬ 
bulent  velocity  profiles  one  would  expect  to  find  in  a  0.635-cm  diameter  pipe,  with  sea  water  at 
25°  C^,  flowing  at  60.26  cm/sec  (Re  =  4060).  Note,  as  a  result  of  increased  turbulent  mixing,  the 
mean  velocity  levels  in  turbulent  flow  change  much  less  (profile  is  flatter)  than  laminar  flow 

^  Discussion  and  photographs  of  the  stabilizing  effect  of  a  contracting  inlet  can  be  found  in  Rouse  (1938). 

®  It  has  been  empirically  observed  that  pipe  flow  remains  laminar  for  Reynolds  numbers  less  than  about  2000 
(White,  1974)  no  matter  how  rough  and  noisy  the  entrance  conditions  are  made.  Laminar  pipe  flow,  after  taking 
great  care  in  minimizing  the  inlet  and  background  disturbance  levels,  has  been  maintained  until  Reynolds  numbers 
of  around  100,000  (Pfenninger,  1961). 

^  Temperature  throughout  the  experiments  reported  herein  ranged  between  15°C  to  25  °C.  This  difference  in 
temperature  results  in  a  20%  change  (decrease  with  increasing  temperature)  in  kinematic  viscosity,  v. 

Consequently,  whatever  the  temperature  chosen  for  a  particular  figure,  the  resulting  pipe  flow  curves  are  fairly 
representative  of  the  entire  experimental  flow  range  tested  herein. 
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across  the  central  core  of  the  pipe.  The  wall  shear  stress  supporting  this  turbulent  velocity  profile 
is  about  two  and  one  half  times  the  laminar  velocity  profile  for  the  same  Reynolds  number. 

As  in  the  laminar  case,  a  relation  between  the  nondimensional  pressure  drop  and  mean  flow 
(i.e.,  X  vs  Re)  can  be  found.  This  empirical  relation,  established  by  Blasius  (1913)  is 

I  =  0.316/Rei/4  (2.10) 

or  equivalently,  via  equations  (2.1b),  (2.5),  and  (2.10), 

Xwall  turb  =  0.04pUavg2/Rel/4  a  (2.11) 

It  can  be  shown  (Schlicting,  1979)  that  equation  (2.11)  is  related  to  the  velocity  profile  specified 
in  equation  (2.9).  In  figure  2.4a  the  theoretical  laminar  and  empirically  derived  turbulent  rela¬ 
tions  for  X  v^'  Re  are  compared  to  tap  water  data  obtained  in  various  pipes  of  different  lengths 
and  diameters.  As  the  ratio  of  entrance  length  (Lg)  to  pipe  diameter  (D)  increases,  either  by 
increasing  Lg  or  decreasing  D,  the  high  laminar  Re  flow  becomes  more  fully  developed.  Note 
also  the  excellent  agreement  between  tap  (figure  2.4a)  and  ocean  (figure  4)  water  data. 

Although  the  dimensionless  variable  X  decreases  with  increasing  Reynolds  number,  it  is 
important  to  realize  that  both  dp/dx  and  twaii  increase  with  flow  speed.  This  is  illustrated  in  fig¬ 
ure  2.4b  where  values  of  wall  shear  stress  and  average  flow  speed  for  fully  developed  flow  in  a 
0.635-cm-diameter  pipe,  with  seawater  at  25°C  flowing  through  it,  are  presented.  Wall  shear 
stress  values  increase  monotonically  with  flow  speed  in  both  laminar  and  turbulent  flow,  exhibit¬ 
ing  a  discontinuous  increase  at  transition.  As  seen  in  figure  2.4b,  when  transition  occurs  at  an 
average  pipe  velocity  of  89.3  cm/sec  (Re  ~  6000)  an  accompanying  increase  in  the  total  shear 
stress  of  almost  3-1/2  fold  is  found  throughout  the  pipe.  With  transition  occurring  at  higher 
Reynolds  numbers,  the  accompanying  increase  in  shear  stress  becomes  larger. 

2.4  LENGTH  SCALES  OF  TURBULENT  EDDIES 

In  fully  developed  turbulent  pipe  flow  the  range  of  eddy  length  scales  can  be  inferred.  The 
largest  turbulent  eddy  length  scale  is  of  the  order  of  the  pipe  radius.  The  smallest  turbulent  eddy 
length  scale,  where  inertial  and  viscous  forces  are  nearly  equal,  is  of  the  order  on  the  Kolmorgo- 
rov  length  scale,  Lk-  The  Kolmorgorov  length  scale  is  obtained  from  dimensional  analysis  after 
making  the  reasonable  assumption  that  the  dissipation  rate  per  unit  mass,  e,  and  the  kinematic 
viscosity,  v  (v  =  p/p),  govern  the  small  scale  motion  (Tennekes  &  Lumley,  1972).  The  Kolmogo¬ 
rov  length  scale,  Lk,  is  defined  as 

Lk  =  (v3/e)1/4  (2.12) 

Experiments  measuring  the  energy  spectra  of  grid  turbulence  (Van  Atta  &  Chen,  1969;  Jayesh  & 
Warhaft,  1992)  have  shown  a  cut-off  close  to  the  wavenumber  corresponding  to  Lk. 

In  fully  developed  pipe  flow,  the  average  dissipation  is  equal  to  the  work  performed  by  the 
wall  shear  stress  acting  on  the  fluid  surface  over  a  displacement  corresponding  to  the  average 
velocity  of  the  flow  (Bakhmeteff,  1936).  Consequently  for  turbulent  flow  between 
Re  =  2300-20000,  where  equation  (2.10)  is  valid,  the  average  dissipation,  £avg,  can  be  esti¬ 
mated  (Davies,  1972)  as 
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Eavg  =  0.16(Re)-0-25(Uavg)3/D. 


(2.13) 


After  substituting  equation  (2.13)  into  equation  (2.12),  one  arrives  at  an  expression  (Rohr  et  al., 
1990)  for  Lk, 

Lk  =  1.58D(Re)-0-69  (2.14) 

Near  the  pipe  wall  a  better  approximation  for  e  can  be  made.  Because  the  velocity  is  zero  every¬ 
where  on  the  pipe  boundary  (no  slip  condition),  at  the  pipe  wall  only  derivatives  with  respect  to 
the  pipe  radius  are  not  equal  to  zero.  Consequently 


e(r-R)  ^  v[(du/dr]2 

(2.15a) 

and  via  equation  (2.7)  with  u  v 

-»0 

E(r^R)  ~  ’t^wall^/(M'P)- 

(2.15b) 

Therefore 

LK(r-^R)  ~  v/(Twal]/p)^^^- 

(2.16) 

Figure  2.5  contains  various  estimates  of  the  Kolmorgorov  length  scale  over  the  flow  range  of 
interest  (pipe  i.d.  =  0.635  cm  and  temperature  =  20.6°C).  The  solid  and  long  dashed  lines  are 
calculated  from  equations  (2.14)  and  (2.16)  which  respectively  employ  average  and  wall  dissipa¬ 
tion  values.  The  short  dashed  line  is  from  a  derivation  by  Davies  (1972),  who  uses  a  dissipation 
representative  of  eddies  near  the  pipe  axis.  As  can  be  seen  in  figure  2.5,  regardless  of  the  esti¬ 
mate  used  for  e,  Lk  is  almost  always  less  than  the  size  of  bioluminescent  dinoflagellates,  the 
smallest  of  which  are  typically  about  30  pm  in  diameter  (Morin,  1983).  When  the  turbulent 
eddies  are  less  than  the  size  of  the  organisms,  they  should  be  much  more  effective  in  distorting, 
as  opposed  to  convecting,  the  organisms. 

Two  issues  concerning  the  above  development  and  interpretation  of  Lk  should  be  discussed. 
First,  the  dissipation,  e,  is  not  uniform  across  the  pipe,  but  decreases  rapidly  from  the  wall 
(Hinze,  1975;  Laufer,  1954).  But  because  e  enters  into  the  equation  for  Lk  as  e^^^,  an  approxi¬ 
mate  value,  such  as  Egyg,  E(r-»0)5  or  E(r_R)  may  be  sufficient.  Second,  significant  shear  energy 
may  only  exist  in  turbulent  scales  10  times  and  greater  than  Lk  (Lazier  &  Mann,  1989).  Accord¬ 
ingly,  Lk  would  have  to  be  about  3  pm  to  ensure  adequate  stimulation  for  the  smallest  biolumi¬ 
nescent  dinoflagellates  present.  As  seen  in  figure  2.5,  higher  Reynolds  number  flows  are  neces¬ 
sary  to  achieve  this. 

2.5  CAVEATS 

It  has  been  tacitly  assumed  throughout  this  analysis  that  the  presence  of  the  organisms  do  not 
significantly  effect  the  pipe  flow  field  described,  and  remain  roughly  homogeneously  distributed 
throughout  it.  While  we  believe  these  assumptions  are  reasonably  justifiable,  there  are  several 
flow-organism  interactions  worth  noting: 

1.  Gyrotactic  focussing.  Kessler  (1985)  has  reported  that  for  small  laminar  flow  rates 
(Re  — 10,  Uavg  —  0.1  cm/sec  for  D  =  0.635-cm)  the  swimming  direction  of  algal  cells  can  be 
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directed  by  compensating  gravitational  and  viscous  torques.  In  downward  pipe  flow  an  initially 
uniformly  distributed  cell  population  becomes  focused  along  the  pipe  centerline. 

2.  Inertial  migration.  Suspensions  of  neutrally  buoyant  spheres  under  laminar  flow  condi¬ 
tions  (Re  =  3-30)  have  been  found  to  migrate,  due  to  the  inertia  of  the  fluid,  to  an  equilibrium 
position  between  the  pipe  wall  and  centerline  (Segre  &  Silberg,  1962a,  b).  This  equilibrium  posi¬ 
tion  has  been  predicted  (Schonberg  &  Hinch,  1989)  to  move  toward  the  wall  with  increasing 
Reynolds  number. 

3.  Turbophoresis.  Particle  drift  velocities  arise  in  turbulent  flow  when  the  statistical  proper¬ 
ties  of  the  velocity  fluctuations  vary  in  space,  resulting  in  a  migration  of  particles  from  regions 
of  high  fluid  turbulence  to  regions  of  low  fluid  turbulence  (Young  &  Hanratty,  1991). 

4.  Drag  and  lift  forces.  Particles  close  to  the  wall  experience  an  augmented  response  to  fluid 
turbulence  (Rizk  &  Elghobashi,  1984). 

5.  Effect  of  crossing  trajectories.  The  trajectories  of  particles  heavier  than  the  suspending 
medium  will  tend  to  cross  the  trajectories  of  fluid  elements,  resulting  in  a  reduced  dispersion 
(Wells  &  Stock,  1983). 

Assuming  the  bioluminescent  organisms  follow  the  fluid  flow,  thereby  ignoring  the  afore¬ 
mentioned  complexities,  is  believed  to  be  a  reasonable  first  approximation  for  the  following 
reasons: 


1.  For  the  most  part,  the  pipe  flows  pertinent  to  this  investigation  are  several  orders  of  mag¬ 
nitude  greater  than  the  swimming  speed^  of  the  organisms. 

2.  The  ratio  of  organism  to  fluid  density  is  significantly  smaller  than  the  particle  experi¬ 
ments  which  exhibited  anomalous  dispersion^. 

3.  Concentrations  of  the  organisms  are  believed  to  be  much  less  than  necessary  to  effect 
flow.  Direct  numerical  simulations  of  particle  dispersion  in  decaying  isotropic  grid  turbulence 
have  shown  (personal  communication.  Said  Elgobashi)  that  10“^  volume  fractions  must  be 
exceeded  for  the  particles  to  significantly  effect  the  flow  field.  This  is  much  greater  than  the  con¬ 
centrations  found  in  the  samples  which  were  inspected. 

4.  Finally,  throughout  the  present  experiments  the  friction  factor  vs  Reynolds  number  data 
for  tap  and  seawater  are  (except  for  transition)  essentially  identical,  suggesting  that  the  flow 
fields  are  similar. 


Even  if  one  knew  everything  about  the  flow  field  in  which  bioluminescence  occurs,  there  is 
still  the  inherent  ignorance  of  the  organisms’  prestimulus  history.  This  is  unavoidable  regardless 

^  In  general,  single  celled  organisms,  typically  between  2  and  1000  pm  long,  can  sustain  swimming  speeds  up  to 
100  body  lengths  per  second  (Yates,  1986).  For  the  organisms  of  interest  {Gonyaulax polyedra  is  about  45  pm 
long  and  41  pm  wide  [Morita  &  Johnson,  1974]),  this  relationship  would  predict  swimming  speeds  less  than  1 
mm/sec,  in  agreement  with  measurements  by  Kamykowski  et  al.  (1992). 

^  When  Young  and  Hanratty  (1991)  measured  the  turbulent  motion  of  100-pm  diameter  glass  spheres  (density  = 
2.42  g/cm^)  suspended  in  water  flowing  through  a  pipe,  they  found  excellent  agreement  with  Laufer’s  (1954) 
classic  pipe-flow  measurements  made  without  particles.  The  density  of  six  marine  dinoflagellates  studied  by 
Kamykowski  et  al.  (1992)  were  between  1.06  and  1.09  g/cm^.  (Note,  some  dinoflagellates  can  achieve  positive 
buoyancy  in  seawater.  Noctiluca  have  a  large  vacuole  in  which  they  concentrate  ammonia  and  exclude  potasium 
and  divalent  anions  to  reduce  density  (Vogel,1981).) 
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of  the  flow  field  chosen  (see  for  example  Latz’s  (1994)  comments  on  the  prestimulus  associated 
with  start  up  in  Couette  flow).  To  test  the  effect  on  bioluminescence  intensity  and  color  spectra 
due  to  transit  in  a  hose,  Losee  and  Lapota  (1981)  pumped  San  Diego  Bay  water  at  0.38  1/s 
through  different  lengths  of  2.54-cm  inside  diameter  (i.d.)  hose^®.  A  viewing  chamber  fitted  with 
quartz  windows  was  positioned  at  the  end  of  the  hose  and  configured  to  provide  additional  turbu¬ 
lent  stimulus.  The  average  measured  bioluminescent  intensity  increased  by  only  22%  when  the 
hose  length  was  decreased  from  120  to  6  m.  The  bioluminescent  color  spectra  showed  little 
change.  It  was  concluded  that  the  longer  hose  caused  little  additional  bioluminescent  degrada¬ 
tion.  Later  measurements  (Lieberman  et  al.,  1987)  showed  that  further  shortening  of  the  hose, 
from  6  to  1  m,  resulted  in  an  additional  increase  in  bioluminescence  of  up  to  30%.  It  is  antici¬ 
pated,  based  on  the  present  investigation,  that  an  even  greater  (possibly  as  much  as  a  factor  of 
two  to  three)  increase  in  bioluminescence  would  be  observed  at  the  hose  inlet.  The  initial  stimu¬ 
lation  associated  with  the  pipe  inlet  is  impossible  to  avoid.  However,  it  appears  that  in  the  region 
of  fully-developed  turbulent  flow  studied  here,  where  the  flow  field  stimuli  is  well  characterized, 
some  quasi-steady  state  evolves  in  which  the  bioluminescent  response  is  not  strongly  dependant 
on  downstream  position.  It  is  in  this  region  that  the  present  study  is  principally  concerned.  It  is 
assumed  that  what  is  learned  here  about  the  fundamental  nature  of  the  response  of  luminous 
organisms  to  hydrodynamic  stimuli  will  be  relevant  to  all  flows  in  general. 


Because  of  the  high  Reynolds  number  and  the  abruptness  of  the  entrance,  it  is  reasonable  to  assume  (via  Eq.  2) 
that  fully-developed  turbulent  flow  was  established  within  the  first  few  meters  of  hose. 
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(R-r)/R 


Figure  2. 1 .  Experimentally  observed  (  ■  ,  □  ;  Laufer,  1954)  distribution  of  the  normalized  turbulent  shear  stress, 
across  a  pipe  in  fully  developed  flow.  Turbulent  shear  stress  levels  have  been  normalized  by  the  wall  shear 
stress,  allowing  different  Reynolds  number  experiments  to  be  compared.  The  solid  line  shows  the  distribution 
of  the  theoretical  values  of  the  total  shear  stress  across  the  pipe.  Except  near  the  pipe  wall  (r  ^  R,  pu'v'  ^  0;  not 
shown)  the  total  shear  stress  is  essentially  entirely  turbulent  in  nature. 
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(b)  Re  =  10500,  beginning  of 
transitional  flow. 


(c)  Re  =  3200,  laminar  flow 


(a)  Re  =  1 3500,  beginning  of 
turbulent  flow. 


Figure  2.2.  Photographs  of  dye  streaks  taken  in  the  0.635-cm  i.d.  pipe 


Laminar 

Turbulent  (Blasius  1  91  3) 

Patel  and  Head  (1  969,  Re=201  5) 
Patel  and  Head  (1  969,  Re=4060) 


Figure  2.3a.  Normalized  velocity  profiles  in  fully  developed  pipe  flow.  Data  from  Patel  and  Head  (1969)  for 
Re  =  2015  (laminar)  and  Re  =  4060  (turbulent).  Solid  line  is  the  theoretical  laminar  velocity  profile.  Dashed  line  is 
the  empirically  derived  (1/7  power  law;  Blasius,  1913)  velocity  profile. 


It 


sec 


r  (cm) 


Figure  2.3b.  Laminar  (solid  line)  and  turbulent  (dashed  line)  velocity  profiles  in  fully  developed  pipe  (0.635-cm 
i.d.)  flow  with  seawater  at  25°C,  for  Re  =  4060. 
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Reynolds  Number 

Figure  2.4a.  Representative  fresh  water  friction  factor  Reynolds  number  data.  As  expected  for  high 
Reynolds  number  laminar  flow  (equation  2a),  with  increasing  hJY)  the  discrepancy  between  measure¬ 
ments  and  theory  (solid  line)  decreases.  Comparison  with  similar  seawater  data  in  figure  4  is  excellent. 
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Wall  Shear  Stress  (dynes/ cm^) 


Reynolds  Number 

Figure  2.5.  Estimated  Kolmogorov  length  scales  (L^)  as  a  function  of  Reynolds  number  for  turbulent  flow  of  2rC 
seawater  through  a  0.635-cm  i.d.  pipe.  The  difference  between  curves  result  from  using  different  dissipation  values, 
i.e.,  average  dissipation  (solid  line;  equation  2.14),  dissipation  at  the  wall  (long  dash;  equation  2.16),  and  dissipa¬ 
tion  throughout  the  central  core  of  the  pipe  (short  dash;  Davies,  1972).  A  cutoff  for  the  size  of  the  smallest  biolumi- 
nescent  organisms  anticipated  to  be  present  is  included  for  comparison  (dashed  line). 
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APPENDIX  3:  INTERPRETING  PMT  MEASUREMENTS  IN  FULLY 
DEVELOPED  PIPE  FLOW  -  THE  ROLE  OF  THE  MEAN  VELOCITY 


An  increase  in  PMT  levels  with  greater  flow  through  the  pipe  may  result  from  both  increas¬ 
ing  hydrodynamic  stimulus  and/or  increasing  mean  velocity,  Uavg-  The  explicit  PMT  dependence 
on  mean  velocity  refers  here  to  the  convective  effect  of  the  flow  on  the  flux  of  luminescent 
organisms,  and  excludes  the  flow’s  stimulatory  aspects.  Individual  flash  measurements  of  biolu- 
minescent  dinoflagellates  have  been  represented  as  where  lo  is  the  peak  intensity 
expressed  as  counts  per  second,  and  y  is  the  decay  rate  of  the  flash  (Seliger  et  al.,  1969;  Widder 
et  al.,  1993).  Modeled  as  such,  the  direct  contribution  of  the  mean  velocity  towards  increasing 
PMT  levels  will  depend  on  the  ratio  of  the  organism’s  residency  time,  L/Uavg  (where  L  is  the 
length  of  pipe  within  the  PMT’s  field  of  view),  to  the  organism’s  flash  decay  rate. 

Bathyphotometers.  -  The  relationship  between  an  organism’s  flash  kinetics  and  its  PMT  resi¬ 
dency  time  has  been  previously  derived  for  bathyphotometers  (Widder  et  al.,  1993)  and  are 
included  here  for  the  sake  of  comparison.  Each  flash  in  the  bathyphotometer  is  assumed  to  be 
initiated  at  the  beginning  of  the  PMT’s  aperture.  The  time  integrated  flash  fraction,  Ip  (counts), 
measured  by  the  PMT  of  a  single  organism  is  the  following; 

r  L/Uavg 

Ip  =  I  lQe~*/‘^dt  =  IqY(1  -  (3-1) 

Jo 

If  (L/Uavg)/Y  «  1,  then  equation  (3.1)  reduces  to 

Ip  =  Io(L/U..g).  (3.2) 

The  mean  intensity,  Imean^  is  proportional^  to  the  product  of  the  average  signal  each  organism 
contributes  and  the  number  of  organisms  passing  through  the  PMT’s  aperture  per  second.  Thus 
for  (L/Uavg)/Y  «  T 


Ijnean  ^  Io(L/Uavg)r'UavgA  ^  ulgLA,  (3-3) 

where  n  is  the  concentration  of  luminous  organisms,  and  A  is  the  pipe  area.  No  explicit  depen¬ 
dence  between  Imean  and  Uavg  appears  in  equation  (3.3).  This  results  because  while  doubling 
Uavg  doubles  the  flux  of  potential  flashes,  the  duration  of  each  flash  within  the  PMT  aperture  is 
reduced  by  one  half.  An  implicit  dependence  remains  however,  since  the  probability  of  stimulat¬ 
ing  an  organism  (and  possibly  the  intensity  and  duration  of  its  flash)  will  depend  on  the  level  of 
flow  agitation  which  is  a  function  of  Uavg- 

If  (L/Uavg)/Y  »  equation  (3.1)  reduces  to 

~  (3-4) 

^  Not  necessarily  equal  because  not  every  organism  passing  the  PMT  is  excited. 
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and  Imean  becomes 


Imean  <x  iQynUavgA. 


(3.5) 


There  is  now  an  explicit  linear  dependence  between  Imean  and  Uavg-  This  results  because  each 
flash  under  the  present  conditions  is  recorded  essentially  in  its  entirety. 

Fully  developed  pipe  flow.  -  Where  the  flow  is  fully  developed,  hydrodynamic  stimulation 
remains  constant  (in  a  statistical  sense)  along  the  flow  direction.  Hence,  flash  initiation  is  equally 
likely  to  occur  anywhere  along  the  pipe.  Note,  hydrodynamic  stimulation  is  not  the  same  across 
the  pipe,  as  larger  shear  stresses  are  found  in  the  slower  flow  near  the  pipe  walls.  Nevertheless,  it 
is  convenient  to  choose  Uavg  ^s  a  representative  velocity  of  the  stimulated  organism,  and  discuss 
the  consequences  of  this  assumption  later. 

The  delay  time,  x,  is  introduced,  so  that  the  product  (Uavg)(x)  determines  the  location  along 
the  pipe  where  the  flash  begins.  The  integrated  flash  fraction  of  organisms  stimulated  within 
(Ipi)  and  upstream  (If2)  of  the  PMT’s  aperture  are  respectively 


■Wu 


'2/  '^avg 


^F1  ” 


loC-b-^Vvdt  =  IoY(1  -  e^/Ye(-L2/U-g)/Y),  (3.5) 


and 


iF2 


r  L2/Ujvg 

’ 

■'U/u.vg 


IqC  =  iQye^A^gC  Li/Uavg)/Y  —  q(  (3.7) 


The  average  intensity  per  flash,  Ipavg,  is  calculated  by  integrating  x  from  0  to  L2/Uavg,  and  then 
dividing  by  L2/Uavg.  Thus, 


•  L,/U, 


^Favg  [ 


H" 


r  L2/Uavg 

•'L./Uavg 


lFidx]/(L2/Uavg)  (3-8) 


which  after  integration  becomes 

Ipavg  =  YIoKLj  -  Li)/L2  -  Y(>Javg/L2)(e-<‘-./n...)/T  -  eC-La/U-J/V)].  (3.9) 

If  (L2/Uavg)/y  and  (Li/Uavg)/Y  «  1,  then 

Ipavg  =  Io(Li  -  Lj)/(2U.vgL2),  (3.10) 

and 

Imean  ^  Iq  (L2-Lj)/(2UavgL2)nUavgA  oc  nAIo(LH2y(2L2).  (3.11) 

Note  that  equation  (3.11)  contains  no  explicit  dependence  on  the  convective  velocity.  Compari¬ 
son  with  the  corresponding  (Li  =  0,  L2  =  L,  and  (L/Uavg)/Y  «  1 )  bathyphotometer  results 


3-2 


shows  that  Imean  defined  by  equation  (3.11)  for  fully  developed  pipe  flow,  is  1/2  of  Imean  defined 
by  equation  (3.3)  for  the  bathyphotometer.  The  factor  of  1/2  arises  because  the  organisms  in  fully 
developed  pipe  flow  are  equally  likely  to  flash  anywhere  within  the  PMT  aperture. 

If  (L2/Uavg)/Y  and  (Li/Uavg)/Y  »  1,  then 

If..6  “  YIo(L2  -  Li)/L2>  (3-12) 


and 


imean  Yio((^2  (3.13) 

An  explicit  linear  dependence  between  PMT  levels  and  the  convective  velocity  now  appears. 
Comparing  again  to  the  corresponding  (i.e.,  (L/Uavg)/Y  »  1)  bathyphotometer  results,  it  is  found 
that  equation  (3.13)  reduces  to  equation  (3.5).  Here  the  agreement  is  exact  because  in  both  cases 
essentially  the  full  flash  is  recorded  regardless  of  where  it  is  initiated  within  the  PMT’s  field  of 
view. 

For  the  0.635-cm  i.d.  pipe  flow  experiments  described  earlier,  and  assuming  Li  =  65  cm, 

L2  =  70  cm  and  y  =  0.1  s,  the  PMT  dependence  on  the  convective  velocity  is  found  (figure  3.1) 
to  be  nearly  linear  throughout  the  entire  flow  regime  of  interest  (Re  <  20,000).  Note  under  these 
conditions,  Li/Uavg/Y  »  1  and  L2/Uavg/Y  »  1-  In  regards  to  using  Uavg  as  the  average  convec¬ 
tive  velocity  of  stimulated  organisms,  note  that  a  smaller  value  which  would  be  more  appropriate 
for  the  highly  stimulatory  flow  occurring  near  the  wall,  only  makes  the  approximation, 
(Li,2/Uavg)/Y  »  better. 
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L]  =  65  cm,  L2  =  70  cm,  v=  0.0104  cm^/sec  (25.6°C).  Decay  rate  of  flash  (y)  is  assumed  to  be  0.1  sec.  Limiting  behavior  for 
(L,  2/Uavg)/Y«  1  (equation  3.1 1,  solid  line)  and  (L,  2/Uavg)/Y»  1  (equation  3.13,  dashed  line)  are  included. 
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